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We present new methods for the evaluation of one-loop tensor integrals which have 
been used in the calculation of the complete electroweak one-loop corrections to e + e~ — > 
4fermions. The described methods for 3-point and 4-point integrals are, in particular, 
applicable in the case where the conventional Passarino-Veltman reduction breaks down 
owing to the appearance of Gram determinants in the denominator. One method consists 
of different variants for expanding tensor coefficients about limits of vanishing Gram 
determinants or other kinematical determinants, thereby reducing all tensor coefficients 
to the usual scalar integrals. In a second method a specific tensor coefficient with a 
logarithmic integrand is evaluated numerically, and the remaining coefficients as well 
as the standard scalar integral are algebraically derived from this coefficient. For 5- 
point tensor integrals, we give explicit formulas that reduce the corresponding tensor 
coefficients to coefficients of 4-point integrals with tensor rank reduced by one. Similar 
formulas are provided for 6-point functions, and the generalization to functions with more 
internal propagators is straightforward. All the presented methods are also applicable if 
infrared (soft or collinear) divergences are treated in dimensional regularization or if mass 
parameters (for unstable particles) become complex. 
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1 Introduction 



Future high-energy colliders, such as the LHC and the ILC, will allow us to search for 
new physics and to test the Standard Model of the electroweak and strong interaction with 
high precision. Various interesting processes naturally involve many particles in the final 
state, where "many" means three, four, or more particles. Such processes often proceed 
via one or more resonances that subsequently decay, or they represent an irreducible 
background to such resonance processes. In order to exhaust the potential of future 
colliders, precise theoretical predictions including strong and electroweak corrections to 
many-particle processes are mandatory. 

The calculation of radiative corrections to complicated processes poses a number of 
problems. Besides the huge amount of algebra, the appearance of unstable particles, and 
the integration of the mult i- dimensional phase space, a numerically stable evaluation of 
the loop integrals is an important ingredient. In this paper we are concerned with the 
calculation of one-loop integrals, including those with five and six external legs. The 
generalization from six to more external legs is straightforward. 

Pioneering work in the calculation of one-loop integrals was performed by Veltman 
and collaborators. Together with 't Hooft, he provided compact explicit expressions for 
the basic one-loop integrals, the scalar 1-point, 2-point, 3-point, and 4-point integrals [UJ, 
which have been completed later by other authors Elaborating on an idea of Brown 
and Feynman together with Passarino he provided systematic formulas that allow to 
reduce all tensor integrals with up to four internal propagators to the basic scalar integrals 
[Hj. These methods are basically sufficient for the calculation of radiative corrections to 
processes with four external particles for non-exceptional configurations. Nevertheless, 
in the sequel some improvements and additions have been worked out. Van Oldenborgh 
and Vermaseren constructed a different tensor basis that allows to concentrate some of 
the numerical instabilities into a number of determinants [Ej. Ezawa et. al performed 
the reduction using an orthonormal tensor basis [E]. A reduction in Feynman-parameter 
space, which is equivalent to the Passarino- Veltman scheme, is used in the GRACE 
package [Ej. 

The main drawback of the Passarino- Veltman reduction and variants thereof is the 
appearance of Gram determinants in the denominator, which spoil the numerical stability 
if they become small. In processes with up to four external particles this happens usually 
only near the edge of phase space, e.g. for forward scattering or on thresholds. For 
the special cases where a Gram determinant is identically zero, alternative reduction 
procedures have been devised by Stuart and collaborators [ |HJ E| (see also Ref. [ITU]). 
However, in processes with more than four external particles, Gram determinants also 
vanish within phase space, and methods for the calculation of tensor integrals are needed 
where Gram determinants are small but not exactly zero. In Ref. [ [H] such a method has 
been devised by constructing combinations of iV-point and (N — l)-point scalar integrals 
that are finite in the limit of vanishing Gram determinants and using this limit if the 
Gram determinant becomes small. 

On the other hand, alternative tensor reduction schemes have been developed using 
different sets of master integrals. Davydychev could relate the coefficients of one-loop 
tensor integrals to scalar integrals in a different number of space-time dimensions [ [J^] , 
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and Tarasov found recursion relations between these integrals [EI]- These methods have 
been further elaborated by different groups [ [TU E3 EH H3 E3 . In this approach all 
one-loop tensor integrals can be reduced to finite 4-point integrals in (D + 2) dimensions 
and divergent 3-point and 2-point integrals in D dimensions. Numerical instabilities in 
this reduction, which are also due to small Gram or other kinematical determinants, have 
been investigated in Ref. [EI] for the massless case, and a systematic improvement by 
an iteration technique has been proposed. While numerically stable analytic expression 
for the basic integrals are available for the massless case, these turn out to be hard to 
construct for the massive case. Therefore, one typically reduces these basic integrals to 
the usual scalar integrals or, in particular for vanishing Gram determinants, calculates 
them by numerical integration [I19j. 

Other algorithms, which are based on recursion relations similar to Passarino-Veltman 
reduction and applicable irrespective of the number of external legs, have been presented 
in Refs. [ BH . These algorithms do not completely avoid the appearance of inverse 
Gram determinants. 

It was realized already in the sixties by Melrose that scalar integrals with more than 
four lines in the loop, i.e., 5-point and higher-point scalar integrals, can be reduced to 
scalar integrals with less internal propagators in four dimensions [1221- These methods 
were subsequently extended and improved by several authors [ [TTJ [T21 EH EH EH 1221 
I2H 1213 I2E1 I2Z] and generalized to dimensional regularization in Refs. [ H3J |2H1 EH] ■ In 
Ref. [ 26J, a method for the reduction of 5-point integrals that completely avoids inverse 
leading Gram determinants has been worked out. Recently, a similar reduction has been 
found that even reduces 5-point tensor integrals to 4-point integrals with rank reduced 
by one [I19j. In all these approaches 5- and higher-point tensor integrals are reduced to 
tensor integrals with less internal propagators. 

Various approaches have been proposed that use numerical integration of loop integrals 
and are, thus, complementary to most of the methods mentioned so far. In the approach 
of Ref. [ I3l)j . which has been elaborated for general one-loop integrals with up to six 
external legs, the Feynman-parameter integrals are rewritten in such a way that they can 
be numerically integrated in a stable way. A fully numerical approach to calculate loop 
integrals by contour integration was proposed in Ref. [Ell- A semi- numerical approach 
that relies on the subtraction of UV and infrared divergences has been advocated in 
Ref. [E2|- A different semi-numerical method makes use of the fact that all tensor one- 
loop integrals can be expressed in terms of one- and two-dimensional parameter integrals 
which are suitable for numerical integration [ I33j . A numerical method based on multi- 
dimensional contour deformation has been proposed in Ref. [E3- Finally, Feynman- 
parameter integrals have been numerically performed with a small but finite "ie" from 
the propagator denominators and a subsequent extrapolation e — > in Ref. [ I34| . So 
far, none of these methods has proven their performance in calculations of higher-order 
corrections for processes with more than four external particles. In practice, one can 
still expect problems with the numerical stability of the algebraic reduction to standard 
forms in specific regions of phase space and with the speed of the underlying numerical 
integration of the basic loop integrals. 

In this paper we describe methods that have actually been used in the calculation of the 
electroweak corrections to e + e~ — > 4fermions [EH], i-e., in the first established calculation 
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of the complete one-loop electroweak corrections to a process with six external particles. 1 
In this approach, 6-point integrals are directly expressed in terms of six 5-point functions, 
and the 5-point integrals are written in terms of five 4-point functions. While we used the 
methods described in Refs. [(221121] and Ref. [|2H] m the original calculation in this 
paper we describe improved methods for the reduction of 6-point and 5-point integrals 
which have meanwhile been implemented in the code for the electroweak corrections to 
e + e + — > 4 fermions and which further improve its performance in numerical stability and 
CPU time. The 3-point and 4-point tensor integrals are algebraically reduced to the 
(standard) scalar 1-point, 2-point, 3-point, and 4-point functions as described below. For 
1-point and 2-point integrals explicit numerically stable results are used. 

In more detail, the 3-point and 4-point functions are reduced to scalar integrals ac- 
cording to the Passarino-Veltman algorithm if no small Gram determinants appear. This 
is the case for most points in parameter space. If a small Gram determinant appears, the 
reduction of 4-point to 3-point or 3-point to 2-point functions is done differently. Here we 
have worked out two alternative calculational methods (referred to as "rescue systems" 
in Ref. [ESI)- O ne method makes use of suitable expansions of the tensor coefficients 
about the limit of vanishing Gram determinants. This is achieved in an iterative way 
and requires to calculate (N — l)-point functions of higher degree compared to the usual 
Passarino-Veltman reduction. 2 Finally, again all tensor coefficients can be expressed in 
terms of the standard scalar 1-point, 2-point, 3-point, and 4-point functions. In practice, 
we use the first two to three terms in the expansions and we have to introduce different 
expansions for different regions of parameter space. In the second, alternative method we 
evaluate a specific tensor coefficient, the integrand of which is logarithmic in Feynman 
parametrization, by numerical integration. Then the remaining coefficients as well as the 
standard scalar integral are algebraically derived from this coefficient. This reduction 
again involves no inverse Gram determinants; instead inverse modified Cayley determi- 
nants appear. In this approach, the set of master integrals is not given by the standard 
scalar integrals anymore. For some specific 3-point integrals, where the modified Cayley 
determinant vanishes exactly, analytical results have been worked out that allow for a 
stable numerical evaluation. 

The paper is organized as follows. We summarize our conventions and useful defini- 
tions in Section |21 The evaluation of 1-point and 2-point tensor integrals is summarized 
in Sections El and 01 respectively. In Section El we provide several methods for the re- 
duction of 3-point and 4-point tensor coefficients and describe their actual application 
to e + e~ — > 4/ in Section 15.71 In Section 15.81 we consider UV and infrared divergences 
in detail and conclude that the proposed methods are valid independent of method for 
infrared regularization. The reduction of 5-point and 6-point tensor integrals to integrals 
with smaller rank and smaller number of propagators is detailed in Sections El and re- 
spectively. In App. |X1 we list the UV-divergent parts of one-loop integrals that enter the 
reduction formulas. Appendix El describes a treatment of singular 3-point integrals based 

1 The GRACE-loop collaboration has recently reported on progress towards one- loop calculations for 
2 — > 4 particle processes. Using the methods described in Refs. [0123, first results on e + e~ — * !/i>HH 
have been shown at conferences and a status report on e + e _ — > /i^D^ud has been given in Ref. [I37j. 

2 A similar idea, where tensor coefficients are iteratively determined from higher rank tensors has been 
described in Ref. [H8|. 
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on analytical methods. Finally, we discuss alternative reductions of 5- and 6-point tensor 
integrals in Apps. O and |Dl respectively. 



2 Conventions and notation 

One- loop tensor iV-point integrals have the general form 

rpN,n\...iJ,p {ph m _ u m0; _ _ _ ; m7v _ i} = ^tl f d D q ^■■■ q " P (2.1) 

m z J N Ni . . . iyjv-i 

with the denominator factors 

N k = (q + p k ) 2 -m 2 k + ie, k = 0, . . . , N - 1, p = 0, (2.2) 

where ie (e > 0) is an infinitesimally small imaginary part. For P = 0, i.e., no integration 
momenta in the numerator of the loop integral, (|2.1j) defines the scalar iV-point integral 
T^. Following the notation of Ref. [H] we set T 1 = A, T 2 = B, T 3 = C, T 4 = D, T 5 = E, 
T 6 = F. Throughout we use the conventions of Refs. [I2UI2E] to decompose the tensor 
integrals into Lorentz-covariant structures. 

In order to be able to write down the tensor decompositions in a concise way we use 
a notation (similar to the one of Ref. [H]) in which curly braces denote symmetrization 
with respect to Lorentz indices in such a way that all non-equivalent permutations of the 
Lorentz indices on metric tensors g and a generic momentum p contribute with weight 
one and that in covariants with n p momenta p^ (j = 1, . . . ,n p ) only one representative 
out of the n p \ permutations of the indices ij is kept. Thus, we have for example 

{P- ■ •/'};' V:/,'.'" = Pi! ■ --Pip, 

{9P\i: P = 9^+9^+9^1, 

{9pp)TJ: = <T<K 2 + g^plPl + g^P-A + 9 vp pIpI + g pa pX + g^fiA 
{ggY vpa = 9 pu g pa + g up g» a + 9 Pil 9 va - (2.3) 

This definition is unique up to the selection of the representative permutations of the mo- 
menta. For our calculation this does not matter, since the covariants are always contracted 
with quantities that are totally symmetric in the indices ij. In fact in our calculation the 
definition is equivalent to a normalization of the sum of the n p \ covariants with a factor 
1/rip!; in this case the third line of (|2.3J) would contain 12 instead of 6 terms on the r.h.s. 
We decompose the general tensor integral into Lorentz-covariant structures as 

[§] N-l 

1 —2.^1^ \9_-£P- ■ ■PJi 2n +i...ip 1 a.^i 2 „+i...ip 

n=0 i 2n+ i,...,i P =l n ^ 
N-l N-l 

= E p£---pt p Tl. ip + E {gp---p}tt p T?oi 3 ...i P 

ii,...,ip=l i3,...,ip=l 
N-l 

+ E {99P ■ ■ ■ i'Y^X:" r mm,:,...ir + ■■■ 

is,...,ip=l 
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+ < 



( N-l 

J2{9---9P}i 1 p - 4lP Tl. 0ip , for P odd, 



{g...g}K~»rTl. , for P even, 



(2.4) 



where [-P/2] is the largest natural number smaller or equal to P/2. For each metric tensor 
in the Lorentz covariant the corresponding coefficient carries an index pair "00" and for 
each momentum pi r it carries the corresponding index i r . 

For tensor integrals up to rank five the decompositions more explicitly read 

AT-1 N-l 



ii=l 



U,«2=l 



rpN,llvp 
rpN,/j,upa 

rj-\N,flUpUT 



N-l N-l 

E WA + EIspK, 

21,12,13=1 11=1 

AT-1 N-l 

E ^J&O^** + E {9VV}tZ a T^ il2 + {^}^T lo, 

il,«2,«3,*4 = l 11,12=1 

AT-1 AT-1 



00iii2«3 



*1,*2,*3,*4,*5 = 1 
AT-1 



U,«2,i3=l 



+ e {sMir^ 



OOOOii • 



(2.5) 



ii=i 



Because of the symmetry of the tensor Tj^ p all coefficients T^ Ap are symmetric under 
permutation of all indices. For convenience we assume this symmetry also for indices that 
are zero. 

When reducing a tensor integral T^ + i p , one encounters tensor integrals that are 
obtained by omitting the fcth denominator N^; we denote such integrals „ P (A;). In 
the decomposition of Tj^ (k), k = 1, . . . , N, shifted indices appear which we denote as 



(ir)* 



i r for k > i r , 
i r — 1 for k < i r . 



(2.6) 



After cancelling the denominator N the resulting tensor integrals are not in the standard 
form but can be expressed in terms of standard integrals by shifting the integration 
momentum. We choose to perform the shift q — > q — pi, so that the following iV-point 
integrals appear: 



T^-^(0) = (2?F ^ ( 9 4 ^ 



17T Z 



D q ^ ... q »p 



N ± ■ ■ ■ N N ' 



(q + Pk - Pi) 2 - m 2 k + ie, 



k = l,...,N. 



(2.7) 



Note that the scalar integral TJf = T and the tensor coefficients Tq , T q 00 , ... are 
invariant under this shift. The other coefficients of (k) can be recursively obtained 

as 



nN 

v 0^0 i l2n+l---jp 

In 



(0) 



'iN 

2n 



^(0), i 2n+ i, ...,i P >l, 



T, 



N 

0...0 li2n+2---«P 

2n 



N 

(0) = -T^o, 2n+2 ..., p (0) - £^o™ 2n+2 ...* P (0), i 2n+2l . . . ,i P > 0. 



(2.8) 



The recursion is solved by 



k N-l 



Tp..£i...ii 2n+k+1 ...i p (0) — ( ^) k ^2[j) ^..0ii...ii,i 2n+k+1 -l,...,i P -l(^)j 

1=0 V / u,...,ii=l ^7 



2n fe 



«2n+fe+l) ■ ■ ■ ,ip > 1- 



(2.9) 



We also use the notation 5^ = 1 — 5^, i.e., X^ijO • •) — Hi^ji- ■ an d employ the 
caret to indicate indices that are omitted, i.e., 



rpN — rpN 

ii...i T ...ip h...i r -iir+i~-ip' 

In the reduction formulas for the (N + l)-point functions the Gram matrix 

' 2pipi . . . 2pip N N 

y 2p N pi . . . 2p N p N ) 

appears. Its determinant, the Gram determinant, is denoted by 

A^=detzW, 

and its inverse can be written as 



(2.10) 



(2.11) 



1 ?W 



(2.12) 



(2.13) 



where z\f^ is the adjoint of z\f\ which can be calculated as 



2pip! ... 2p 1 p j _ 1 2p 1 p j+l ... 2p x p N 



(N) 



'-lY+3 



2pi-iPi ■ 
2p i+ iPi ■ 

2p N Pi ■ 



2p i _ 1 p j _ 1 2p i _ 1 p j+1 
2p i+1 p j _ 1 2p i+l p j+l 

2p N p j _ 1 2p N p j+1 



2pi-lPN 

2p i+1 p N 
2p N p N 



(2.14) 



i.e., from a reduced determinant of where the ith row and the jth column have been 
discarded. 
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We introduce a generalization of the adjoint by 
v+j +k+i ggn ^ _ k j ggn ^ _ ^ 

2p 1 p j _ 1 2p 1 p j+1 ... 2pxVi-\ 2pip J+ i 



Z { t N) 



(2.15) 



x 



2pfc+iPi 
2p N pt 



2-Pi-iPj-i 2p i _ 1 p j+1 
2p i+1 pj^i 2p i+1 p j+1 

2p fc _ip 3 -_i 2p fc _ip j+1 
2 Pfc+iJ5j-i 2p fc+ ip i+ i 

2p N pi_ x 2p N p j+1 



2pi_lP/_i 2pi_iJ9; +1 
2p i+ ip/_i 2p i+ ip i+ i 

2p fe _ip/_i 2p fc _i^ + i 

2p fc+ ip;_i 2p k+1 pi +1 

2p N pi_i 2p N p i+1 



2p x p N 

2pi-\PN 

2p i+1 p N 

2p k -iPN 
2p k+1 p N 

2p N p N 



i.e., it is defined from a reduced determinant of where the zth and Zcth rows and the 
jth and Zth columns have been discarded. Moreover, it is defined to vanish for i = k or 
j = I. For the case N = 2, it is given by 



K2) 



&Ofci — OinOkl- 



'(ik)(jl) — °il°kj "ij 



(2.16) 



Expanding the determinant of along the fcth row or the Zth column, respectively, 
it can be written as 



N N 
A (AO _ V- - sr 



(2.17) 



m=l m=l 

where k and / are not summed. Expanding the determinant in (|2.14|) with the Zth column 



replaced by (Z\hi 



Z^k) along the Zth column yields the relation 

N _, 

^ y(N) 7 (N) _ y(N), y(N), 



and analogously 



m=l 



V 



m=l 



(2.18) 



(2.19) 



These imply the equations 



^il ^kj ~ ^ij ^kl 



/AW (2.20) 



and 



N 



E z 



(N) 



r(N) _ y(N) 



m,n=l 



(im)(jn) ran ^i] 

Y~ 7 iXi _ aWa a 7W7W 



2V 



(2.21) 



m,n=l 
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An important special case of the last relation is 



N 



\{N) _ 7 (N) y(N) v- y{N) 7 {N) 7 
A - L lk ^lk + 2^ ^(ImMkn^mk A 



AN) 
J ln ■ 



(2.22) 



m,n=l 



The relations ()2.13|1 - ()2.22|1 are valid for any (not necessarily symmetric) N x N matrix 
Z( N ) with determinant A^ N \ 

We further introduce the (N + 1) x (N + 1) matrix 



(2ml h ... f N \ 
fx 2p x p\ ■ ■ ■ ZpiPN 

\ f N 2p N px . . . 2p N p N ) 



(2.23) 



with 



r 2 2,2 

fk = Pk ~ m k + m 0, 

Its determinant is given by 



k = l,...,N. 



(2.24) 



N 



det(XW) = 2m 2 A^ - £ f n f n Z™ 



n,m=l 



Yqo Yqi ■ ■ ■ Yqn 
Y\q Yu . . . Y 1N 



Y m Y m ... Yj 



NN 



det(y), (2.25) 



i,j = 0,...,N. 



(2.26) 



where 

Yij = ml + m)- (pi -pj) 2 , 

The matrix Y = {Yij) is sometimes called modified Cayley matrix and its determinant the 
modified Cayley determinant [123- Its elements are related to those of the Gram matrix 
via 

^ = Z\f ] -f t - fj + 2ml Y oi = Y i0 = + 2m 2 , i,j = l,..., N. (2.27) 

The vanishing of det(X^) is a necessary condition for the appearance of leading Landau 
singularities [EHf. The adjoint of x\f\ i,j = 0,...,N can be expressed as 



vW _ a(JV) 
^oo — ^ > 



N 



n=l 



N 



n,m=l 



iV. 



(2.28) 



For later use we also consider the generalized adjoint of X^ N \ The relevant part of it is 
given by 



(JV) 

(0i)(0j) 



X 



(AT) 



4J , i,j = l,...,JV, 

AT 



(oi)O'fe) - (ifc)(oi) 



n=l 



i,j,k= 1,...,N. 



(2.29) 



These relations together with ()2.20|) imply 



det(xW)ziP = AWX<P -iff 



det(XW)X (7V) 



Oj > 



3 Evaluation of 1-point functions 

The scalar 1-point integral for an arbitrary complex mass m is given by 



A (mo) = ml 



A + In I ^ ] + 1 



where A is the standard one-loop divergence 

A 2 



4- D 



7e + ln(47r) 



(2.30) 



(3.1) 



(3.2) 



in D space-time dimensions with 7e denoting Euler's constant. The tensor integrals of 
rank 2n (n — 1,2, . . .) are given by 



(3.3) 



where the tensor coefficients are easily evaluated to 



A ...o 



in, 



2n 



2 n (n + 1)! 



(3.4) 



Because of Lorentz invariance obviously all tensors of odd rank vanish. 
4 Evaluation of 2-point functions 

In the following we assume that at least one of the parameters p\ , mo, mi is different 
from zero; otherwise the 2-point integrals identically vanish in dimensional regularization, 



£...(0,0,0) = 0, 

where the dots stand for any Lorentz index or any index of a tensor coefficient. 
Up to rank 3 the 2-point tensor integrals are decomposed as 



(4.1) 



B» = p»B x 



B>» =p^ 1 B 11 +g^B l 



oo ■ 



= rfplrfBm + {gpyrBooi. (4-2) 



The tensor coefficients can be algebraically reduced to scalar 1- and 2-point integrals, 
Aq and Bq, with the Passarino-Veltman algorithm [ H] as more generally described in 
the next section. The corresponding results for tensors up to rank 3 are, e.g., given in 
the appendix of Ref. [I2S]- The algebraic reduction for the coefficients B 00i3ii , which 
correspond to covariants involving the metric tensor, 



Boo 
Boot 



1 

6 

1 r 



A o (0) + hB x + 2m 2 B + ml + ml- \p]\ , 
-A (0) + hB n + 2mlB x - \{2m\ + Am\ - p\) 



etc., 



(4.3) 
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are numerically well behaved. However, the reduction formulas for the coefficients 
corresponding to the covariant p^ 1 ■ ■ -pi P involve a factor l/p\ in each reduction step, 
so that these reduction formulas become numerically unstable for small p\. Owing to 
the simplicity of 2-point integrals it is, however, possible to derive closed expressions for 
these coefficients that are numerically stable for all values of p\. Such a derivation is de- 
scribed below. Assuming the knowledge of the coefficients Bi 1, the remaining coefficients 
-S0...01...1 can be obtained from the recurrence relations 



2n + 2 P-2n-2 



1 



2{P - 2n - i; 

P-2 



A0...0 (0) + /iB ...o + 2p 2 P ...o 



2n P-2n-l 



2n P-2n-l 



2n P — 2n. 



71= 0, 



(4.4) 



or 



-B0...0 1...1 



2n+2 P-2n-2 



2(P + 1) 



^^a^o 1...1 (0) + 2mlBo_o + /1-B0...0 1...1 

2n P-2n-2 2n P-2n-2 2n P-2n-l 

P-2" 



- 2(£> - 4)5^0 

2n+2 P-2n-2 



n = 0, . . . , 



The coefficients v4 ...oi...i(0) are given by 

A0...0 (0) = (-l) p - 2 - 1 i ... (0), 



(4.5) 



(4.6) 



2n P-2n-l 



where t4o...o(0) can be obtained from ()3.4|) . The finite polynomial quantities (P — 4)I?oo... 
can easily be derived by exploiting ()4.5|) for the UV-singular parts; explicit results for 
tensors up to rank 5 are summarized in App. [X] 

We derive the expressions for Px...i by explicitly solving the Feynman-parameter inte- 
gral 

Bi^i = dx {-x) n {A + In /i 2 - In -p\x{l - x) + mjj(l - a;) + m 2 x - ie } . (4.7) 

n 

In the following result we support complex mass parameters; more precisely, the real parts 
of m 2 must be non-negative, the imaginary parts negative or zero. The final results are 
conveniently written as 



Bus = ^ { A + In ( 4) - E /„(**)] 

^ 71+1 [ V m 0/ ,5 J 

with Xfc denoting the solutions of the quadratic equation 

= — p\x{l — x) + rra^l — x) + m 2 x — ie. 
For p 2 = one of the x^ is formally 00. The auxiliary functions 

f n (x) = (n + l)£dtt n \n (l-~) 



(4.t 



(4.9) 



(4.10) 
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can be evaluated in a numerically stable way by choosing one of the two representations 



- h ( 1 -i) + ,E«Tr (4U) 

The first form is numerically stable for intermediate values of \x\ 7^ 0. For x — > 0, f n (x) 
develops a true logarithmic singularity; for x — *■ 1 the logarithm ln(l — 1/x) is suppressed 
because of its prefactor. The second equality in (j4.11j) yields numerically stable results 
for large \x\. In practice, we take the first form for \x\ < 10 and the second otherwise. 
The case where one of the Xk is zero corresponds to mo = and can be easily obtained 
via taking the limit mo — > 0, 

BM, mi ) = t£ { A + ln ( -^W) + 4t " /. (l " ^) ) 
v -v^ / n+1 [ \m( — pf — le J n + 1 \ pf J) 

(4.12) 

For = m\ this further simplifies to 



Bu K, 0, rm) = ^ { A + In f 4) + ( 4 - 13 ) 
n + 1 I V mf y n+1 1 



In the vicinity of the last two special cases one of the Xk becomes small, so that the leading 
(logarithmic) term in f n (xk) cancels against the explicit logarithm in (J4.8|) . Although 
this somewhat worsens the precision of the evaluation, we did not find problems with this 
approach in practice. Nevertheless we have additionally implemented a more sophisticated 
representation of with more branches where such cancellations are avoided. 

In the above derivation we essentially followed the approach described in the appendix 
of Ref. [Hj; the results given there are, however, not applicable to the general case of 
complex masses. 

5 Reduction of 3-point and 4-point functions 

The tensor decompositions of 3-point tensor integrals up to rank 4 and 4-point tensor 
integrals up to rank 5 read explicitly 

c p = E pZA,, c» = E p^Pi 2 c ili2 + <rc M , 

h=l 11,12=1 

c^ p = E PtPiPl z c hi2iz + Eiwircoon, 

ili*2,j3=l ii=l 

c pvpa = E p>i 2 p p i 3 Plc ili2i3ii + E {gppWC c oohi2 + {ggr upa c 000 o, 

U,i2,*3,*4 = l 11,22=1 

(5.1) 
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il=l 



D 



= 1 11,12=1 

3 3 

E PhPi2Pi 3 D hi2i 3 + E {MiTAjOn , 



ii=l 



[11/ pa 



D 



[IVpCTT 



*1>*2 1*3=1 

3 3 
E PhPi 2 Pi 3 Pl D hi2i3U + E {^M^f^Wia + {^j^AxTOO, 

il>*2,*3,*4=l 11,12=1 

3 3 

E PilPi2PhPlPi S D hi2i3i4iB + E {3PPP}iZh D Whi2h 

il>*2,*3,*4,*5=l il,*2,«3=l 

+ E{^P}r ,7r ^oooon- (5-2) 
n=i 

Because of the symmetry of the tensor Tj^^ all coefficients C^.^p, and D^ip are 
symmetric under permutation of all indices. To be specific, in the following we give the 
reduction formulas for the 4-point functions, i.e. N = 4. To obtain the corresponding 
results for 3-point functions one has to perform the substitutions 

CL->S„, D... -> C.„, Z( 3 )^Z( 2 \ A^A^, X( 3 >-»X< 2 >, iV^3, (5.3) 

and similar obvious substitutions. 

5.1 Conventional Passarino— Veltman reduction 

The one-loop tensor integrals can be reduced to scalar integrals recursively by inversion 
of systems of linear equations [ H]. The inhomogeneity of these equations consists of 
coefficients of lower rank. The equations of this system are obtained by contracting 
Tf^ ^ p with the (N — 1) external momenta p^ 1 and for P > 2 also by contraction with 
the metric g^ 1 ^ 2 . Contracting (|2.1j) with p^ 1 and using 



2p k q = N k -N -f k 



(5.4) 



each of the first two terms on the r.h.s. of (J5.4j) cancels exactly one propagator denominator 
of Pf^Tj^ ^ p and the third term is proportional to T^ 2 . Likewise the contraction of 
(12. 1J) with g^ 1 ^ 2 yields a factor q 2 in the numerator of g^^Tjl , which can be written 



as 



iVo + ml 



(5.5) 



The Nq term cancels the first propagator, the second term leads to the tensor . 



This yields 



2p?T» 



^jt 2 ...Rp(^) -^i 2 -Rp(0) /fc^)i2-pr' 

V ' (0)+mg2^ 



fJ.lfJ.2rfN _ rp 

y ± Ulfl2- ■fJ.P t 



■ fi 3 ...fip\ 



- (U3— (Up 



(5.6) 
(5.7) 



Note that for T^ 3 „ p (0) a shift of the integration momentum — > q^—pi has to be done 
in order to achieve the standard form (j2.1j) . The tensor integrals with shifted momenta 
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T / ^~j [p (0) are defined in (|2.7|) . Inserting the Lorentz decompositions (15.2)1 into ()5.6|) and 
(15. 7)) . the desired recurrence relations can be read off by comparing coefficients. 
From (J5.6j) we obtain 



ii...ip 



Ski 2 ...i P = C(i 2 ) k ...(i p ) k (k)5 ki2 . . . 8 kip - Cj 2 ,..j p (0) - f k D 

N ~ 1 (3) P 
m=l r=2 

fc = l,...,iV-l, z 2 ,...,z P = 0,...,iV-l, (5.8) 

and from ()5.7|) 

Sm 3 ...i P = 2C i3 ...i p (0) + 2m 2 Di 3 ,„ ip 

N-l P 

= Z^D nmi3 ... ip +2(D + P-2 + J2Ko)D 00i3 ... lp , 

n,m=l r=3 

i 3 ,...,i P = 0,...,N-l, (5.9) 

where the matrix Z^ 3 > is defined in (12.11)1 . Equations ()5.8)) and ()5.9)1 can be solved for the 
coefficients of D^ 1 -^ as 



-Dnn; 



00l3...lj: 



2(3 + P-iV) 



-2(D - 4)A)0i 3 ...iP + Ci 3 ...^(0) + 2mjj£> 



2, ,. 

13... i P 



(5.10) 



JV-l 

"I" E^ fn-D n i 3 ...i p 
n=l 

N-l / P \ 

A 1 .., P = E^ 3 ')^ C..i,- 2 E^^,.i.., P , ii^O. (5.11) 

n=l \ r=2 / 



The relations ()5.10jl and 1)5.11)1 determine D i± Ap in terms of D il ip _ 1 and 3-point 
functions. Using these relations recursively, all coefficients of 4-point functions can be 
expressed in terms of 3-point functions and the scalar 4-point function D . The finite 
polynomial quantities (D — 4)Dooi 3 ...i P can easily be derived by exploiting ()5.10)1 for the 
UV-singular parts; explicit results for tensors up to rank 7 are summarized in App. [X] As 
explained in Section IB~8*l IR divergences do not occur in D 00i3 ip . More explicit formulas 
for all tensor functions up to rank 5 are given in the appendix of Ref . [ [2H] ■ 

Figure n illustrates the Passarino-Veltman reduction scheme for 4-point integrals in a 
plane of tensor coefficients where the rank of the tensor increases by going down in the 
rows and the number of index pairs "00" increases by going to the right in the columns. 
The steps in the algorithm are indicated by arrows that show which coefficient is deduced 
from previously calculated ones. The numbers close to the arrows correspond to the step 
number which is identical to the rank of the tensor coefficients to be calculated; the labels 
"a" , "6" , etc. give the order in which the coefficients within a step are calculated. 

Equation (j5.11)l becomes numerically unstable if Z^> is nearly singular, i.e., if the 
Gram determinant is close to zero. Reduction schemes for this case are described in 
Sections I5.3H5.6L 
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Figure 1: Schematic illustration of conventional Passarino-Veltman reduction. 

5.2 Alternative Passarino— Veltman-like reduction 

An alternative to the conventional Passarino-Veltman reduction can be obtained as 
follows. Equations ()5.8|) and ()5.10p can be written as 



( 2m 2 


h h h \ 




( A 2 . 


,A 




( A 2 . 


dp 


fi 


2pipi 2pip 2 2pip 3 






..i P 






..i P 


h 


2p 2 Pi 2p 2 p 2 2p 2 P3 




D212 


..i P 


D 2i2 


..i P 


V h 


2p 3 Pi 2p 3 p 2 2p 3 p 3 J 




\D 3i2 


..i P ) 




\Dsi 2 


..i P } 



(2(D + P- N)D mi2 ...ip - C l2 .. Ap (0) \ 



qP 

°li 2 ...ip 

Soi„ ;„ — 2 Z) r=2 02i r D mi2 I 



2i 2 ...ip 

P o v^-P 

3«2...ip 



r=2 "-Uf^OOia... ir» .if 
ip 

r.-.ip 



2 Y, r =2 ^3i r D 00i2 - ir 



12, 



i P = 0,...,N-l, (5.12) 



where on the r.h.s. the matrix A® defined in (|2.23|) appears and the following abbrevia- 
tions are introduced, 



qP 



Ski 2 ...i P + /fcA 2 ...ip- 



C (i2)k-(ip)k( k )^ki 2 ■ --Skip ~ C i2 .. Ap (0) 

Multiplying (J5.12)) with the matrix A® from the left, we obtain 
det{X®)D i2 „. ip = A (3) [2(4 + P - N)D 00i2 .., ip + 2(D - A)D 00i2 ... ip - C, 



(5.13) 



12— lp 



(0) 



(3) 



and 



2V-1 
n=l 



det(X^)D ni2 ... lp = m 2(4 + P - N)D 00i2 ... lp + 2(D - 4)Z> 



^ p — 9 V A D 

u ni 2 ...i P * n V 00i 2 -V- -»P 



r=2 



(5.14) 



00i2-"«p 



.(0) 
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N-l 



n=l 

Equation (|5.14|) yields D, 

2(4 + P-N)A^D 00l2 ... ip 



Sni 2 ...i P 2 E 5ni r D t 



r=2 



H ^ 0. (5.15) 



ooi 2 ..i P in terms of /^ m , , , D h _ ip , and 3-point functions, 



-2A< 3 )(/J - 4)£W.i P + A (3) C l2 ...i P (0) 



AT-1 



+ det(X( 3 ))A 2 ...i P - E^< 



(3) 
On 



n=l 



00i2-. .i r ...ip 



r=2 



(5.16) 



and thereafter (|5.15jl yields D^...^. Using these relations recursively, all coefficients of 
4-point functions can be expressed in terms of 3-point functions and the scalar 4-point 
function D . While the final results are of course identical to those of the usual Passarino- 
Veltman reduction, the order in which the different coefficients are calculated is different. 
As a consequence, this recursion can, in some cases, be numerically more stable than the 
conventional Passarino-Veltman reduction, in particular, if all the quantities A^ 3 \ xj^ , 
and X$ become small. 

For the tensor coefficients up to rank 3 the reduction formulas explicitly read 



N-l 



2(5 - N)A^D 00 = -2A^(D - 4)Ax> + A^O) + det(X^)D - £ (5.17) 



n=l 



det{X^)D h = X$[2(5 - N)D Q0 + 2{D - A)D 00 - C (0)] + E (5.18 



N-l 



(3) Al 



n=l 



2(6-iV)A (3) £) 1 



00«2 



-2A^(D - 4)A) 0l2 + A^C t2 (0) + det(X^)D l2 

N-l 

~ E \Sni 2 — 25 ni2 D 00 , 

n=l 

det(X^)D tll2 = X$[2(6 - N)D 00i2 + 2(D - 4)D 00h - C i2 (0)] 



N-l 



+ E^ 



(3) 



71=1 



S„ in — 2S n i 2 Doo 



m 2 



%iM ^ 0, 



2(7 - A0A (3) £> 



-2A^\D - 4)D 00l2l3 + A< 3 >C iai8 (0) + det(X^)D 

N-l 



E x < 



(3) 
On 



n=l 
'(3)[ 



S 



2d n i 2 Dooi 3 — 25 n ioD, 



ni-j, J -^00i2 5 



det(X( 3 ))A 1J2i3 = *$[2(7 - N)D mi2lA + 2{D - A)D 00i2l:i - C l2iz 

N-l 



(3) 



((V.2'3 



25 nin D, 



ni2 -LsOOis 



2d n i 3 Dooi 2 



(5.19) 



(5.20) 



(5.21) 



«i,i 2 ,i 3 ^ 0. (5.22) 



'2 ''3 



n=l 



Note that (|5.21|) holds also for 12 = 23 = 0. 

The 3-point tensor coefficients that result from omitting iVo in the 4-point integrals 
are defined according to (|2.8|) or more explicitly 



C n (0) = ^-iW, h = 2,...,N-l, 



N-l 



n=2 



c 1 (o) = -ECn(o)-c (o), 
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(5.23) 
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Figure 2: Schematic illustration of alternative Passarino-Veltman reduction. 
C hi2 (0) = 6^-1^-1(0), h,h = 2,...,N-l, 

N-l 

Ciii(O) = - E C««(0) - C n (0), ii = 1, . . . , N - 1, (5.24) 

n=2 

Ctli2i3(0) = Cii-l,i2-l,i3-l(0)) ^2, ^3 = 2, . . . , N — 1, 

Ci hi2 (0) = - E ^W ) - C W2 (0), M2 = 1, • • • ,N - 1. (5.25) 

Figure 121 illustrates the alternative Passarino-Veltman reduction scheme for 4-point 
integrals in the plane of tensor coefficients similarly to Figure of the previous section 
for the conventional variant. 



5.3 Reduction with modified Cay ley determinants 

Equation (|5.12|) can also be exploited directly to calculate tensor coefficients of lower- 
rank from higher-rank tensors. Specifically, the coefficients -D^...^ with %\ 7^ for tensors 
of rank P are expressed in terms of the coefficients Dooi 2 ...i P for tensors of rank (P + 1). 
This means, (|5.12|) recursively expresses tensor coefficients D^ . ip in terms of C functions 
and of a single coefficient -D ...o which results from D i± ip upon replacing all non-zero 
indices i& by "00". For sufficiently high tensor rank P, viz. P > 2N — 4, the integrand 
of the Feynman parameter integral of -D0...0 involves only polynomials and logarithms 
of the integration parameters x/. Such integrals are numerically well behaved, because 
singularities appearing in logarithms can be safely treated numerically. The explicit form 
of the Feynman-parameter integral for the general coefficient with P > 2N — 4 is 
given below. 

In summary, equation (|5.12|) provides a method for deducing all tensor coefficients 
D^ ^ip (including the standard scalar integral Do) from C functions and the numerically 
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evaluated coefficient D of tensor rank 2P. This procedure does not involve the inverse 
of the Gram determinant A^ 3 ), as it is the case in the two versions of Passarino-Veltman 
reduction described in the previous sections. However, the method involves the inverse 
of the modified Cayley determinant det(X( 3 )), so that it becomes unstable if det(X^ 3 ^) 
becomes small. It is also interesting to note that the numerically evaluated coefficient 
-D0...0 enters this reduction with a prefactor A' 3 '. Thus, this method becomes particularly 
precise if A^ 3 ) is small, where Passarino-Veltman reduction is unstable, because the error 
in the numerical calculation of -D0...0 is suppressed in this case. Note, however, that both 
the reduction of this section and Passarino-Veltman reduction become problematic if both 
A® and det(X^) are small. 

For tensor coefficients up to rank 3 the reduction formulas explicitly read 

JV-l 

det(X^)D 0000 = A< 3 )[2(9 - N)D 000000 + 2{D - 4)A)ooooo - C ooo(0)] + £ X$S 5 n0000 , 

n=l 

(5.26) 

JV-l 

det(X( 3 ))A) = A® [2(7 - N)D 0000 + 2(D - 4)A)ooo - Gx>(0)] + E ^S 3 n00 , (5.27) 



n=l 

JV-l 



det(X^)D = A( 3 )[2(5 - N)D W + 2(D - A)D 00 - C (0)] + E X^S*, (5.2 

71=1 



det(V( 3 ))Ax>00H = X®[2(9 - N)D 0000Q0 + 2(D - 4)/J 00 oooo - C 00 oo(0)] + E X hl^ 



N-l 

V-(3) c 

nOOOO' 
n=l 

(5.29) 

N-l 

det(X^)D 00ll = X\%2(7 - N)D 0000 + 2{D - 4)/J 00 oo - Coo(0)] + E ^S^noo, (5-30) 

n=l 

"(3)[ 



det(X^)D 00lli2 = Xgj[2(8 - N)D QQ00l2 + 2{D - 4)D 00Q0i2 - Q 



00i 2 



N-l 



+ E X hn Sn00i 2 ~ 25 n i 2 Z)QO0O , H,^ 7^ 0) (5.31) 



n=l 



Finally, D it , D ili2 , and D ili2i3 are obtained from (|5.18j) . (|5.2(J|) . and (|5.22|) . respectively. 
Thus, all 4-point tensor coefficients up to tensor rank 3 can be recursively deduced from 
-Doooooo and 3-point coefficients. 

Figure HU illustrates the reduction scheme for 4-point integrals up to rank 3 in the 
plane of tensor coefficients similar to the previous sections. The steps of the reduction 
now proceed from right to left, starting with a basis integral -D ...o with as many index 
pairs "00" as the finally aimed tensor rank, i.e., for rank 3 with -D ooooo- in each step we 
get all coefficients of at least one rank lower with one index pair "00" less than in the 
previous steps. 

Generically the Feynman-parameter integral for reads 

N-A-2k 



1 



2 fe (2 + k — N)\ 



N-l \ /N-l 



Tl, = ofc , , 7 m n / dx, 6 l-J2a lXl J2x m j\2+k-N 



X 



j=0 u j \ 1=0 / \m=0 / 

£;>iV-2, (5.32) 



2+k-N -i /A-\f\ f N ~ l 

E V" +2k U> 

n=l 11 \ A 4 / \m=0 
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Figure 3: Schematic illustration of the reduction with modified Cayley determinants. 



with the shorthand 



fN-l 



'N-l 



'N-l 



A = A(x , . . . , Xjv-l) = X lPl ~ X) X m X n(pl ~ 



nr 



(5.33) 



\m=0 



,n=0 



The real parameters a\ appearing in (|5.32|) are widely arbitrary; they only have to fulfil 
the constraints oci > and X^q 1 a i > 0- For the numerical evaluation of the Feynman- 
parameter integral it is convenient to take the uniform choice ot\ = 1, in which case the 
integral runs over the (N — 1) -dimensional unit simplex (T/v-ij 



N-l j-l 

x = l-^2x h 0<xj<l-^2xi, j = l,...,N-l, 
1=1 1=1 



(5.34) 



The integral representation (|5.32|) is valid both for real and complex masses. 

Specifically the integrals for C ooooo an d -Doooooo, which are needed for tensors of rank 3, 
are given by 



a 



""° 000 = 2880 ( A + D [ S ^ 2 + p i + p 2 + s ^{p\ + Pi) + p\p\ 

of 2 , 2 2 , 2 2\ 

- 3(m si2 + p 1 m 2 + p 2 m i) 

- 6[si 2 (m 2 + ml) +pj(ml + m 2 ) + pl(ml + m 2 )] 
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i 1 E r~ 4 I 4 i 4 i 22, 22, 22 

+ 15 [m + m 1 + m 2 + m Q m l + m m 2 + m 1 m 2 



-D oooooo — (A + 1, 



•^0l2 + Sia + S 23 + Pi + ^2 + Ps) + ^( m l + m l + m l + m l) 



1 /• , 3 M-ie 



d ^i n /__| (5.36) 

with the shorthands 

S12 = (pi - P2) 2 , S13 = (pi - p 3 ) 2 , s 23 = (p 2 - P3) 2 - (5.37) 

For an efficient numerical integration of these integrals we use a fortran code based on the 
DCUHRE algorithm [ 39], as included in the CUBA library [ 40j. The UV-divergent 
parts are integrated analytically in order to ensure exact cancellation of the singularities. 

As mentioned above, the procedure described in this section becomes unstable if 
det(X < - Ar_1 ' ) ) becomes small. The basis integrals are still safely calculated via the 
numerical integration, but using the described relations to deduce the remaining coeffi- 
cients accumulates an instability in each step that turns an index pair "00" into a non-zero 
tensor index or that eliminates an index pair "00" . This accumulation of an instability 
can be suppressed by extending the set of basis integrals. For instance, the 3-point tensor 
coefficients can be deduced from the coefficients Coo, Coooo, and Coooooo, which all 

have logarithmic integrands in their Feynman parametrizations, upon using the above re- 
lations only once. 3 If det(X^ N ~~^) is not small, we prefer to deduce all tensor coefficients 
from one basis integral (e.g., -Doooooo f° r D ili2i3 ), because no instabilities accumulate and 
the recursion preserves relations among the tensor coefficients, which are less accurately 
valid if several coefficients are calculated numerically. 

If det(X( N ~^) = 0, the described procedure is not applicable. For instance, this is 
the case for 3-point functions that are either soft or collinear singular. Such cases are 
much simpler than the case with general kinematics, so that they can be treated more 
directly. For processes with light external fermions only, det(X^ -1 ^) is zero only for 
3-point functions (N = 3) where a photon or a gluon is attached to an external fermion. 
A fully analytic treatment of these cases, which admits a numerically stable evaluation, 
is described in App. |Bj this method can be extended to similar cases that appear in other 
processes. 

Finally, we remark that the method of this section is somewhat related to the fully nu- 
merical procedure advocated in Ref. [ EOj . There, a method is described how the Feynman- 
parameter representation of one-loop integrals is, upon partial integrations, transformed 
into integrals with logarithmic integrands, which are then treated numerically. The oc- 
curring algebraic coefficients that express the original integral in terms of logarithmic 
integrals are related to the coefficients of the inverse of the matrix X^ N ' introduced in 
this paper. In fact we have verified that the reduction of the scalar integral Cq to loga- 
rithmic integrals leads to the same results as our equation (|5.28|) for iV = 3 [see (|5.3|) ]. 

3 Note that the Feynman-parameter integral of D Q0 is not logarithmic, so that the calculation of D ili2i3 
from -Dqooo an d -Doooooo requires the use of the recurrence relations twice. 
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Therefore, like in our approach, also in the approach of Ref. [ 30j, the cases with small or 
vanishing modified Cayley determinant det(X^ N ~ 1 ^) require a special treatment. More- 
over, we emphasize that we treat only one basis integral numerically, while the procedure 
of Ref. [ED] in general involves more numerical integrals. 



5.4 Reduction for small Gram determinant 

Let us now derive a reduction scheme that can be used if the Gram determinant A^ 3 ^ 
becomes small, but without changing the set of basis integrals, which are thus still the 

_____ ~s /q\ 

standard scalar integrals A , B , C , D . Multiplying ()5.8|) with indices nix ■ ■ - ip by Z] n 
and summing over n yields 

jv-i / p \ 

X§D h ... ip = - £ Z$ S^ ip - 2 £ 6 nir D 00h .j_ ip + A«%., p (5.38) 

n=l V r=l / 

for arbitrary j = 1, . . . , N — 1 and i r = 0, . . . , N — 1. In order to arrive at this form, 
()2.17|) and ()2.28|) have been used. As long as at least one of the quantities , defined in 
(|2.28jl . is large compared to A^ 3 \ ()5.38|) can be used to determine Di x _ ip from D 00ii j r ■ 
up to terms that are suppressed by the factor A®. 

In order to obtain -D 00il j- i , we consider for arbitrary k, I ^ 



N-l 

A (3) D Hl3 .., P = ' p 



N-l / _ N-l z \ 

= Y \^kl %ij + Y ^(kn)(lm)^nj Z im \ D iji3 „ Ap , (5.39) 
i,j=l \ n,m=l J 

where f!2.21|) has been used. The first term on the r.h.s. can be reduced with ()5.9|) . the 
second term on the r.h.s. upon using ()5.8|) twice. Collecting terms containing Dooii...i P 
and making use of ()2.19|) and (|2.21|) . we obtain 



2 [6 + P - N + £ KoJ Z { $D mi _ Ap = -2(D - A)Z^ D 00il .. Ap - A^D 

N-l 

I V(3) qP+2 i / 7(3) qP+2 7(3) qP+2 \ 

°00ii. ..i P y^nl D nkh...i P n U D nni 1 ...i P ) 

n=l 

N-l ~ r P 

~ Z {kn){lm) fn S m t L ... ip + 2 Y $ni r ^ m 00n..i r ..ip ~~ fnfmP > i 1 .. 



kli\...ii 



n,m=l 

P P 

2^(/n^mir + fmfinir)P , QOi l ...i r ...ip ~ 4 ^ ^ni r ^mi s Dooooi 1 ,,2 r ...t a ...i I 

f=\ r,s=l 



(5.40) 



which holds for arbitrary k, I = 1, . . . , N — 1 and ix, ■ ■ ■ , ip = 0, . . . , N — 1. Together 
with ()5.38|) this equation allows to iteratively determine the tensor coefficients of 4-point 
functions in terms of 3-point functions for small Gram determinant A^. If the 3-point 
functions are known up to rank P, all 4-point tensor coefficients up to this rank can be 
determined recursively up to terms of order A® from these equations by putting all terms 
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involving A^ to zero. Inserting these results back into the r.h.s. of ()5.38|) and (|5.40|) for 
the terms proportional to A( 3 ), all 4-point tensor coefficients up to rank (P — 1) can be 
determined up to terms of order (A^ 3 )) 2 , and so on. Finally, the scalar 4-point function 
is iteratively determined up to terms of order (A^ 3 ^) p+1 . In order to improve numerical 
stability, we can choose j in (|5.38|) such that xff is maximal, and k and I in (|5.40|) such 

that is maximal. For A^ 3 ^ = this reduction scheme essentially corresponds to the 
one proposed in Ref. 

For the lowest tensor coefficients the explicit results read 

N-1 

X®D = -Y,zgs 1 n + A®D j , 



n=l 
K3); 



2(6- N) ZtfDu = -2(D - 4)^ } A)o - A®D kl + Z { k fs oo 



N-1 N-1 _ 

' y^nl °nk "hi °nn) Z^ 



n=l 



71,771=1 



(kn)(lm) 



(3) Q 2 
fnS m — fnfmDo 



A 0j U ii 



N-1 



z ZZ^(S 2 mi -25 nil D 00 )+A^D Jll 



71=1 



2(8- N) zi?D mi = -2(D - 4)Z$D mi - A^D klll + Z$S 3 mi 



(3) 



K3) c3 



N-1 N-1 

El 7(3) A3 _ 7(3) A3 \ _ 7(3) 
l^nZ °rifcii ^fcZ °rmiiJ Z^ Z/ (fcn)(Zm) 
71=1 71,771=1 



f Q 2 _|_ 9A . Q : 
Jn u m i 1 ~ 4u Mi IJ ) 



mOO 



^ (3) n 

^Oj -^1112 



fnfmDi^ 2(f n 5 m i 1 ~\- fm^nii) Dqq 
N-1 



X] ^in ^niii 2 ~~ 2(<5 )1 j 1 Dooi 2 + ^nia-^oon) + A^ D jhi . 2 , 



n=l 



2(8- N) zj^Dwoo = -2(D - 4)zl?D 0000 - A< 3 >£>ooh + ^^oooo 



AT-l 



, / 7VV c4 7(.Jj q4 

"r ^ni °n00fc ^kl D n 



OOn/ 



71=1 



2 (10 - N) zifD 00hi2 = -2(D - 4)Z$D mii2 - A^D kUll2 + Z^S^ 



N-1 _ 

7(3) 

Z^ ^(fcn 
n,m=l 

(3) 



(kn)(lm) 



fnS m QQ f n f m D t 



00 



K3) C 4 



N—l 



f 7P) Q4 A4 
"1 Z^ V^ni D nkhi 2 ^kl °nni 1 i 



(3) A4 



71=1 

AT-l _ 



Z\kn)Um) fn&mixiz + 2($niiS m00 i 2 + 8 n i 2 S mQ0il ) fnfm^iii 



(kn)(lrn) 
n,m=l L 



2(/-n<^mj2 "I - fm&ni2 )-^00ii 
0000 j ^1;*2 



^Oj " ty 21«2«3 



AT-l 

( 3 ) n _ _ V 7 (3) 

Z^ 

71=1 

«2«3" 



•^niiiaM ~ 2(5 n i 1 D 0i 2 i3 + ^7W2 A)Qiii3 + ^ra 3 -DoOii 



+ A®D jhi 



(5.41) 

(5.42) 
(5.43) 



(5.44) 
(5.45) 

(5.46) 



(5.47) 



(5.48 



Figure 0] illustrates a systematic algorithm for this iteration scheme for 4-point in- 
tegrals in the plane of tensor coefficients similar to the previous sections. Thin arrows 
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Step 2 




Step 3 




Figure 4: Schematic illustration of the iteration for small Gram determinants, where thin 
arrows indicate that the relation involves a suppression factor A^ 3 ). In each step the 
boxed coefficients are calculated in the order indicated by the labels "a", "6", etc. The 
nth iteration consists of the following (n + 1) steps: n—>-(n — !)—>... — »1— >0. 



indicate that the relation involves a suppression factor A^ 3 \ At the beginning of the 
iteration all 4-point tensor coefficients as well as the scalar integral Dq are set to zero, 
i.e., no 4-point basis integral is needed. The nth iteration consists of the (n + 1) steps 
n—>-(n — 1)— ► . . . ^1^0 and requires all 3-point coefficients of rank n. Step n with n > 
starts with the coefficient of rank (n + 1) with the highest number of index pairs "00" , i.e. 
with the right-most coefficient in the (n + l)th row in the diagrams in Figure 0] Within 
a step, coefficients for rank (n + 1) are deduced from the right to the left in the diagram; 
only for the last coefficient (which has no index pair "00" ) one has to go one step upwards 
to rank n in addition. After the nth iteration the tensor coefficients -Dj 1 i 2 i3... °f rank n 
without index pairs "00" and all coefficients Dooiii 2 - °f one ran k higher with at least one 
index pair "00" are obtained up to terms that are suppressed by a factor A^ 3 ). Coefficients 
of a rank that is lower by a number m are known up to terms suppressed by [A^]" 1 " 1 " 1 . 
Indicating coefficients that are known up to terms of O (jA( 3 )] m+1 ) with a superscript 
"(m)", the iteration proceeds as follows: 
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Iteration 0: is calculated; all other coefficients are still zero. 

Iteration 1: Step 1 yields Dqq and D^; step yields Dq . 

Iteration 2: Steps 2 to deliver D$ h , DfJ 2 , D$ , , and L>J 2) . 

Iteration 3: Steps 3 to deliver D^, D^ lia , D^ h , D$ , D?' , and 



D ( 



)(3) 

'o 
etc. 



~ (3) 

The reduction method described in this section breaks down if none of the X$ is 
large compared to A^ 3 ) or if all zff become small, since in these cases the iteration does 
not converge. A reduction for small X^j is described in Section 1531 A reduction for 
small Z kl is given in Section l5~El For N = 2 the case of small Z kl is equivalent to small 
zg ) ;foriV = 31 
configurations. 4 



Z^i ; for iV = 3 the case of small covers the case of small apart from exceptional 



5.5 Reduction for small Gram determinant and small modified Cayley de- 
terminant 

If in addition to the Gram determinant A^ 3 ) also all quantities X^ \ j — 1, . . . , N — 1, 
become small, the reduction scheme of Section 15.41 breaks down. As can be seen from 
()2.3()j) . in this case the determinant det(A^ 3 ^) = det(y) of (|2.25|) becomes small, which is 
a necessary condition for the appearance of leading Landau singularities. In this situation, 
we can determine the tensor coefficients as follows. 

For i r 7^ 0, equation ()5.38|) can be rewritten as 

p _ N-l 
2 E ^'V = E Z^S^ Ap + X$D h ... ip - \ :i P f ,,. (5.49) 

r=l n=l 

This allows to determine D 00i ^ f for . . . , ip ^ in terms of 3-point functions as: 

N-l 

ZP^kl Dqq^i = E Z kn $nV..l + X k0 ~ ^D k j_^, 

P-i n= l p P P 

N-l 

2(P - l)Z kl Dqq^jj^ = -^Z kh D 00 j_^ + ^ Z kn S n + L..lh 
P-2 P-l n=1 p_i 

+ xj$D L „ lil -A®D kUlil , h^0,l, 



p-i 



2(P — 2)zlfD 00 L..ii 1 i 2 — 2 Zfcii Dqq i _ j j 2 -2Z k \ > 2 D 00 i,,ij 1 + E Z kl^ + ^ 



N-l 

n I. ..I ii%2 



P-2 P-2 



4 In an alternative approach, one could disregard l|5.40[l and use l|5.38|) also to determine D 00i j - . 

~ (3) . . T 

This reduction method would also work if all Z kl are small. However, in this case, tensor integrals of 
higher rank would be needed. For instance, to calculate Di 1 i 2 i 3 in leading order in A^ 3 ) one would have 
to calculate -Dqooooo and Cqooooo first. 
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1. ..I 1112 

F-2 



A (3) D 



k 1...1 iit2 i 

P-2 



h,i 2 ^ 0,1, 



(5.50) 



and so on, provided that at least one of the Z kl is not small. Again k ^ and / ^ 

~ (3) 

can be chosen such that Z kl is maximal in order to improve the numerical stability. The 
tensor coefficients with more index pairs "00" can be determined by equations that are 
obtained from ()5.50|) by adding additional index pairs "00" to all quantities S and D in 
fl53nj) . 

In order to derive a relation for the calculation of D ix . _ ip we rewrite (15.12)) as 



( fi f2 h \ 

2pipi 2p x p 2 2pip 3 
2p 2 p 1 2p 2 p 2 2p 2 p 3 
V 2p 3 p x 2p 3 p 2 2p 3 p 3 j 

(2(D + 1 + P 

qP+l 

°lh...ip ' 



( -Diji...i P \ 



D 2l 



i—ip 



\ 



qP+l 
°2i 1 ...i P 

^Si\...ip 



\D 3il ...i p j 
2Ef=i 8&..D, 



■ip 



(0) 



2mgA 1 ...i P \ 



OOii ...i r ...ip 



/lAl...i P 



(5.51) 



/ 



After discarding the (j + l)th of these equations, where j = 1, 2, or 3, the remaining three 
equations have the solution 



N-l 



V f 7 (3) n-- 



■ip 



7(3) 



n=l 



2(D + 1 + P - N)D. 



00h...i P 



a 



N-l 



^\in)(jm)f n 



m,n=l 



qP+l 



2^5 mir D, 



00i\...i r ...ip 



i...ip(0) 
fm,Di 1 ...i p 



2uiQDi 1 i p 



r=l 



(5.52) 



Using ()2.28|) this can be written as 

~ (3)n - ■ - zf [2(5 + P - N)D 

N-l 

+ E zg 

m,n=l 



^ij - Ly ii...t.p 



K3) , 

' (* n )(i m ) ^ n 



00ii...ip 

p 



+ 2(L> -A)D 



00ii...ip 



-a 



i\...i P 



(0)] 



S r ni 1 ...i P 2^^5 m i r D, 



00ii...i r ...ip 



r=l 



I ^( 3 ) n 



(5.53) 



which holds for arbitrary i,j 



N — 1 and i±, . . . ,ip 



N — 1. Together 



with ()5.49|) this equation allows to iteratively determine the tensor coefficients of 4-point 
functions in terms of 3-point functions for small Gram determinant A® and small X^f 

~ (3) ~ (3) 

and Xqj as long as at least one of the Xh is not small. Again % and j can be chosen 

~ (3) 

suitably in order to improve the numerical accuracy, e.g. by choosing the maximal X^. 
If the 3-point functions are known up to rank P, all 4-point tensor coefficients up to rank 
(P— 1) can be determined up to terms of order A^ 3 \ xj^f , and X^f from ()5.49|) and ()5.53|) 
by putting all terms proportional to these quantities to zero. Inserting these results back 
into the r.h.s. of these equations, all 4-point tensor coefficients up to rank (P — 3) can be 

] 2 , and so on. Finally, the scalar 



determined up to terms of order [max(|A® |, \, \X^\ 
4-point function is determined up to terms of order [max(|A®|, |A^|, \X { 



: (3) |)][( P + l)/ 2 ]. 
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For the tensor coefficients up to rank 3 the reduction formulas explicitly read 

AT-1 



xWd = Z%> [2(5 - N)D 00 + 2(D - 4) A>o - C (0)] 

AT-1 _ 

+ E Z(in)(jrn)fnS m + X Qj Di, 
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x {3) n 



4^i? ) -Doo«jii 2 
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y {3) d 



n=l 
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(5.54) 
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Figure 5: Schematic illustration of the iteration for small Gram and modified Cayley 
determinants, where thin arrows indicate that the relation involves a suppression factor 
A^ 3 \ -XfcQ , or Xq^ . In each step the boxed coefficients are calculated in the order indicated 
by the labels "a", "6", etc. The nth iteration consists of the following (n + 1) steps: 
7i->(ra- !)-►... ->l-*0. 



jv-i _ 

m,n=l 

+ Dii li2 i 3 . (5.61) 

Figure El illustrates a systematic algorithm for the iteration scheme for 4-point integrals 
in the plane of tensor coefficients similar to the previous sections. Thin arrows indicate 
that the relation involves a suppression factor A^ 3 \ X$, or X$ . At the be ginning of 
the iteration all 4-point tensor coefficients as well as the scalar integral Do are set to 
zero, i.e., no 4-point basis integral is needed. The nth iteration consists of the (n + 1) 
steps n^(n — 1)— » . . . — >1— »0 and requires all 3-point coefficient functions up to rank 
2(n + 1). Step n starts with the two coefficients of ranks (2n + 2) and (2n + 3) that 
have exactly one index pair "00", i.e. which belong to the second column in the respective 
rows, in the diagrams in Figure El Within a step, first the two coefficients are calculated 
that are reached upon omitting the index pair "00" from the starting coefficients; they 
are located in the first column two rows above the starting rows in the diagram. Then 
all coefficients that lie to the right of the starting coefficients are calculated column by 
column. After the nth iteration the tensor coefficients D ili2i3 of ranks 2n and (2n + 1) 
without index pairs "00" and all coefficients D 00ili2 __. of two ranks higher with at least 
one index pair "00" are obtained up to terms that are suppressed by a factor A^ 3 \ X^ \ 
or Xqj . Coefficients of a rank that is lower by a number 2m are known up to terms 
suppressed by [max(|A^|, |A^|, |Xg^|)] m+1 . The iteration proceeds as follows: 



+ ^2 ^{in)(jm)fn 
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• Iteration 0: D$$, D$ h , , and DjJ are calculated; all other coefficients are still 
zero. 

• Iteration 1: Step 1 yields £>S2 x i 2 > £>fii li2 ; 3 , D^} 2 , D|3 ais , D$ 00 , and D$ 00il , step 
yields D&\ D$ k , D«, and D^. 

• etc. 

The reduction described in this section breaks down if none of is large compared 
to A® and , or if all Z^' become small, since in these cases the iteration does 
not converge. A reduction for small zj^ , and thus for small zj^ in non-exceptional 
configurations, is described in Section l5~Tfl If both and all xj^f and xff become 
small, in some cases the alternative Passarino-Veltman reduction of Section I5~2*l works. 



In other cases, none of the discussed reduction methods is really good. However, this 
happens only in exceptional cases, and one of the discussed methods yields at least crude 
results. 5 

5.6 Reduction for small momenta 

(3) 

Finally, we provide a reduction scheme for the case where all Z kl and thus all momenta 
become small. Note that in this case also all of the quantities A^ 3 \ zff, X^, and xj^' 
become small. If the /& are not small as well, we can proceed as follows. We rewrite (|5.8j) 

as 



P N-l 



■ipi 



r=l m=l 

k = l,...,N-l, i 1 ,...,i P = 0,...,N -1, (5.62) 

and ()5.9|) as 

p 

2(4 + P + £ Ko)D 00h ... ip = -2(D - A)D mi ... lp + 2C n .., p (0) + 2mlD h .„ ip (5.63) 

r=l 

N-l 

y , Z nm D nm i 1 ^ i p , i\, . . . , ip 0, . . . , N 1. 

n, m=l 



By using these equations iteratively, we can determine D ix _i p and A)0ii...ip f° r given 3- 

(3) 

point functions for small Z M . If the 3-point functions are known up to rank P, we can 

7 (3) 
J kl ) 

r(3)l 2 



(3) 

determine the coefficients of the 4-point functions with rank P up to terms of order Z k 



those of rank (P — 1) up to terms of order [Z kl '] , . . . , and those of order up to terms 
of order [Z^] p+1 . In order to improve numerical stability, we can choose k such that f k 
is maximal. Note that the structure of ()5.62|) and (|5.63|) is similar to the one of ()5.38|) 



5 An alternative reduction could be derived, by considering X^T^ 1 X^f 1 D i± j ip i , using (|5.53|) and 
inserting l|5.49|) on the r.h.s. of the resulting equation to eliminate D QQi j . . From the obtained 



OOii.-.i, — ip + i ' 

hod wot 

small. However, in this case, tensor integrals of higher rank would be needed. 



relation, all tensor coefficients could be calculated. This reduction method would also work if all Z k ? are 
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and (|5.40|) . In fact, a systematic algorithm for this iteration scheme for 4-point integrals 
is given by Figure EJ if the arrows that point vertically downwards or horizontally to the 
left are omitted. 

Up to tensor rank 3 the explicit formulas read: 



N-l 

A3) 



f k D = Sl-Y: Z^D m , (5.64) 

m=l 

N-l 

8D 00 = -2(D - 4) A)o + 2C (0) + 2m 2 D - £ Z®D nm , (5.65) 

n,m=l 

N-l 

f k D h = S 2 kli - 25 kil D 00 - £ zflD mil , (5.66) 

m=l 

N-l 

12D 00n = -2(D-4)Do 0il +2C r fl (0)+277igA 1 - £ Z nl D nmn, (5.67) 

n,m=l 

N-l 

/fcAii 2 — Skhi 2 ~~ 25fcj 1 Z) oi a ~~ 25fci 2 I?ooii — £ Z km D mili2 , (5.68) 

m=l 

iV-1 

16Aoni 2 = -2p - 4) D 00ili2 + 2C ili2 (0) + 2m 2 D hi2 - £ Z^D nmni2 , (5.69) 

n,m=l 

/fcAii 2 i 3 = Skiii 2 i 3 ~ 25fci 1 DoOi 2 i3 ~~ 25fcj 2 Aoiii 3 — 25fcj 3 Aoiii 2 
JV-1 

~~ E Z krn D mil i 2 i 3 . (5.70) 

771=1 

If also all the /& become small we can rewrite ()5.62|) and ()5.63|) as 

P " iV_1 (3) 

r=l m=l 

fc = l,...,JV-l, n,..., i P = 0,...,JV-l, (5.71) 

and 

p 

2mjjA 1 ...i J , = 2(4 + P + £ 5<ro) Aon...ip + 2(P - 4) Aon... tp 

r=l 

-2C t ,.. lp (0) + £ ^iAmi,...^, i ls ...,ip = 0,...,JV-1. (5.72) 



n,m=l 



By using these equations iteratively, we can determine Ai...i P and A)0ii...tp f° r given 3- 
point functions for small and small The structure of (|5.71j) and (|5.72j) is similar 
to the one of (J5.49|) and (|5.53j) . If the 3-point functions are known up to rank P, we 
can determine the coefficients of the 4-point functions with rank (P — 1) up to terms of 
order max(|Z^|, |/ n |), those of rank (P — 3) up to terms of order [max(|Z^|, |/fc|)] 2 , 
and so on. Finally, the scalar 4-point function is determined up to terms of order 
[max(i4?|,i/ fc |)][(^l. 
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Section 


ED 






EH 


EB 


EH 


method 
type 


PV 
red. 


PV 
red. 


Cayley 
red. 


Gram 
exp. 


Gram/Cay ley 
exp. 


momenta 
exp. 


applicability 




1 r 7\ -in 

\X\JO 


\X\JO 


|Z|->0 
X Oj ^0 
Zm^O 


|Z|->0 

x ti ^ 


Z — >U 


stable for 
\Z |— > ? 


no 


no 


yes 

J 


yes 

J 


yes 


yes 


stable for 
|X|-> ? 


yes 


no 


no 




yes 


yes 


stable for 
? 


no 


no 


no 


no 


no 


yes 


fast ? 


yes 


yes 


no 


yes 


yes 


yes 



Table 1: Summary of features of the reduction schemes for 3- and 4-point tensor integrals. 
The type of the method is either "reduction (red.)" or "expansion (exp.)"; Gram and 
Cayley determinants are generically indicated by \Z\ and respectively. 

5.7 Summary of reduction schemes and application to e+e - — > 4f at one 
loop 

Table ^ briefly summarizes some of the features of the described reduction schemes for 
3- and 4-point tensor integrals. The type of the method, "reduction (red.)" or "expansion 
(exp.)" is indicated in the third row. In the fourth row we summarize the conditions for 
the applicability of the schemes. Conditions that depend on indices i,j,k,l have to be 
fulfilled for at least one choice of these indices. The "yes" and "no" in the last rows 
indicate whether a method is stable or unstable in the corresponding limits or if the 
method is fast in terms of CPU time. A blank entry means that the method can be stable 
or unstable in the considered limit. 

The reduction schemes described above have been successfully applied in the calcula- 
tion of the complete one-loop corrections to the charged-current processes e + e~ — > 4/ as 
presented in Ref. [EH]- As described there, actually two independent calculations of the 
corrections have been carried out employing two different procedures (called "rescue sys- 
tems" there) for the evaluation of the one-loop tensor integrals in the numerically delicate 
kinematical configurations. Both procedures make use of the conventional Passarino- 
Veltman reduction (see Section |5~T|) as long as internal consistency checks prove this 
method to be reliable. If this is not the case, the procedures differ: 

(i) Procedure 1: reduction with modified Cayley determinants and further exception han- 
dling 

If conventional Passarino-Veltman reduction seems not to be trustworthy, since consis- 
tency relations among the tensor coefficients are valid only to very few digits or even vio- 
lated, the method with modified Cayley determinants is used as described in Section 15.31 
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Because of the vanishing modified Cayley determinant this is not possible for the IR- 
singular (i.e. soft or collinear divergent) 3-point functions. Therefore, these cases are 
evaluated as described in App. |H1 yielding perfectly stable results. 

The described procedure fails if both the Gram and the modified Cayley determi- 
nants are very small. In practice, this happens only at a small fraction of events that 
hardly contribute to the e + e~ — > 4/ cross section. However, if this limitation of the 
procedure becomes serious in other cases, the double limit of small Gram and modified 
Cayley determinants can be covered using the method of Sections 15.51 etc.. as it is done 
in Procedure 2. 

(ii) Procedure 2: expansions for small Gram determinants etc. 

The second procedure is based on the two versions of Passarino-Veltman reduction of 
Sections 15. II and l5~2l and on the expansion methods described in Sections l5.4H5.6l If the 
Passarino-Veltman reduction fails, at least the Gram determinant of the corresponding 
integral (or an integral related to a subgraph) is small. The question which of the different 
expansions is most appropriate is decided by estimating the number of valid digits in each 
of the expansion variants; the variant promising the highest precision is taken. 

For the application to the processes e + e~ — > 4/, it turned out to be sufficient to 
implement the expansions for small Gram determinant ( Section I5.4j) . for small Gram and 
modified Cayley determinants ( Section I5.5J1 . and for small momenta ( Section 15.6)1 up to 
tensor rank 4 for 4-point functions and the corresponding formulas for 3-point functions 
up to tensor rank 5. The implementation of the modified procedure for small fk was not 
required. We also did not yet implement the schemes mentioned in footnotes |U and El 

Note that in the one-loop diagrams for e + e~ — > 4/ 3- and 4-point functions appear 
only up to rank 3, i.e. the implemented reductions go beyond taking pure limits of vanish- 
ing determinants. For these processes, the exceptional cases where none of the expansions 
is good appeared only for a very small fraction of events and did not yield sizeable con- 
tributions to integrated physical quantities. 

5.8 UV and IR divergences in dimensional regularization and terms of order 
(D-4) 

In the preceding equations we have kept all terms of order (D — 4) that multiply 
one-loop coefficient integrals. These terms give rise to finite terms in dimensional regu- 
larization if these integrals are divergent in four dimensions. It is convenient to discuss 
UV divergences, which formally result from loop momenta q tending to infinity, and IR di- 
vergences, which arise from finite loop momenta but specific kinematical configurations, 
separately: 

(i) UV divergences 

UV divergences are universal in the sense that the divergent terms in an integral are 
regular functions of the external momenta pk and internal masses m^, but these terms 
do not change if these kinematical quantities approach exceptional configurations (zero 
limits, on-shell configurations, etc.). At one loop, UV divergences generally have the 
form 1/(D — 4) times a polynomial in pk and ntf.. Therefore, the terms (D — 4)T^ ip 
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contained in the above formulas are finite polynomials in p^ and m^. We have listed these 
(D — A)T t N i terms for 1-point functions of arbitrary rank, 2-point functions up to rank 5, 
3-point functions up to rank 7, 4-point functions up to rank 7, and five point functions 
up to rank 6 in App. [XJ 

(ii) IR divergences 

IR divergences at one loop originate either from soft or collinear configurations of a 
loop momentum [ . These type of divergences have the property that they do not 
show up in tensor coefficients with at least one index pair "00", i.e. all tensor structures 
containing at least one factor of the metric tensor are IR finite. 

These fact can be seen by inspecting the Feynman-parameter integrals of the tensor 
coefficients or by the following arguments. A soft singularity results from the limit of 
zero-momentum transfer of a massless particle (q — > 0) between two on-shell particles. 
Assuming that the massless particle correspond to the propagator denominator N , power 
counting in q shows that soft divergences can appear only in the scalar integral, but not 
in tensor integrals, because loop momenta in the numerator render the limit q — > in the 
integral non-singular. Thus, in the general case, where the massless particle corresponds 
to any propagator denominator jV^, soft divergent parts of tensor integrals are always 
proportional to powers of the momenta p^, as can be seen by performing a shift q — > q—pk, 
which maps to N^. A collinear singularity results from the range where the loop 
momentum q is parallel to the momentum pk of a light external on-shell particle that splits 
into two light particles. If a tensor g m . . . q^ p is present in the loop integral, the divergence 
can only show up in covariants that are built up in the singular region. Thus, collinear 
divergences of tensor integrals appear in covariants containing only the momentum p^. 

In the reduction formulas given above the factor (D — 4) appears only in front of tensor 
coefficients A containing at least one index pair "00", which have been shown to 
be IR finite. Therefore, all the reduction formulas are valid without modification if IR 
singularities are regularized dimensionally. All terms (D — 4)Tjj^ 3 i can be taken from 
App. [XJ if more of these terms are needed, they can be easily derived from the reduction 
formulas themselves. 

6 Reduction of 5-point integrals 

In four space-time dimensions, 5-point integrals can be reduced to 4-point integrals. In 
Ref . [ G21 we have given relations that express 5-point tensor integrals of rank P by 4-point 
tensor integrals of rank P (see also App. [0]). This method follows the strategy proposed 
in Ref. [122] for the reduction of scalar integrals and was actually used in the calculation 
of one-loop corrections to e + e~ — > 4/ [EH]. Here we derive formulas that directly reduce 
5-point tensor integrals of rank P to 4-point tensor integrals of rank (P — l). While similar 
results have been presented in Ref. [H2], our derivation is more transparent. 
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We start by considering the determinant 
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(6.1) 

In the first manipulation, we have split the determinant in the second column, and in the 
second we have added the second row of the first determinant to its first row and we have 
moved q a from the first row to the second column in the second determinant. Moreover, 



we have used the definitions ()2.2j) and (|2.24|) . 

In four dimensions, the determinant S vanishes, as can be seen from its defining form, 
because q is linearly dependent on the four momenta pi, i — 1, . . . , 4. Since we want to 
derive a relation that also holds in dimensional regularization we do not use this fact, 
but translate the integral over 8 into a form that has a factor of 0(D — 4) rendering the 
whole contribution zero for finite integrals. Inserting the first form of £ in (j6.1j) into the 
integrand of the tensor integral E^^^ p results in 
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This form can be written more compactly by introducing the four-dimensional metric 
tensor 
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(6.3) 



leading to the result 

JS = 2ml^{g»-g {A f a )E a ^ 
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(6.4) 



The second term is obtained by expanding the second determinant in (J6.2|) along the first 
two rows and the first two columns according to (J2.17j) and using (J2.2U)) . 
Alternatively integrating over the last form of £ in (|6.1|) . we obtain 
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The last determinant can be written as 
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The first equality in ()6.6|) can be easily checked by expanding along the second column. In 
order to explain the second equality, we introduce the Lorentz-covariant decompositions 
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The operation "(g)" isolates all tensor structures in which the first Lorentz index appears 
at a metric tensor; the remaining part of the tensor furnishes the "(p)" contribution. The 
last decomposition in ()6.7|) becomes obvious after performing a shift g — > q — p% in the 
integral. From ()6.7|) it follows immediately that the terms in the second line of (j6.6|) 
that involve [D afll ^^ p (i)]^ p ' , i = 1, . . . ,4 drop out when expanding the determinant along 
the second column, because the resulting determinants vanish. Similarly, the contribu- 
tion proportional to [D afllflp (0) + piaD^ ... MP (0)]^ vanishes after summation over all 
contributions. The remaining terms involving [Dy 3 ' are collected in the quantity 



pam~w(i) = [D a ^~» p (i) - D^ 1 -^ (0)}^ , i = 1,...,4. 
Inserting ()6.6|) into ()6.5|) and expanding the determinants we find 
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where -X"a2v mi * s gi ven m Q2.29J1 . Setting this equal to (|6.4)L we obtain 
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+ 2m 2 A^(g^-g {4 f a )E^-^ 

+ 2 £ X3H [pZ(g a p - (4W ) - V nA€ - 9(0] Erf"*"**. (6.10) 

71=1 

In this result, all inverse Gram determinants have been absorbed in the four- dimensional 
metric tensor, which appears only in the difference (g — gm)- In four dimensions, all these 
terms vanish identically. In dimensional regularization they contribute only if 
involves singularities, i.e., only the singular terms in i£ £ w il -w are relevant. As explained 
in Section IB~Hl IR singularities of fi^ 1 "* appear only in contributions that are propor- 
tional to a momentum p%. These contributions vanish exactly in (j6.10|) as long as the 
external momenta have only non- vanishing components in the four- dimensional subspace. 
UV singularities appear only if P > 4. Therefore, we can omit the last two terms in 
(I6.10|) for P < 4. 6 For P > 4 the inverse Gram determinant that is implicitly contained 
in gu\ can always be cancelled by a prefactor A^. In the last-but-one term of (j6.10|) this 
pref actor is already explicit; for the last contribution it is straightforward to check 7 that 
this factor always arises after symmetrizing the r.h.s. of (|6.1U|) w.r.t. the indices /i, //i, 
. . ., \i P . 

The next step consists in the insertion of the decompositions of tensor integrals into 
Lorentz covariants. Here and in the following we omit the terms involving (g — g^) if 
P < 4. The general tensor decompositions up to rank 5 explicitly read 
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In four dimensions, the covariants involving metric tensors are redundant in these decom- 
positions, since the metric tensor could be replaced by (|6.3j) . By keeping these coefficients 
we can avoid the appearance of explicit inverse Gram determinants in the reduction for- 
mulas. 



6 This result is in agreement with the observation made in Ref. [EH] that in the absence of UV diver- 
gences reduction formulas valid in 4 dimensions remain valid in D dimensions up to terms of 0{D — 4), 
independent of the possible occurrence of IR singularities. 

Contributions to E a/3fJll '" fJ ' p involving p£ vanish, and those involving g^gPM cancel after symmetriz- 
ing w.r.t. the indices /x, hi, . . ., /ip. In terms involving g a @ the surviving (g — guj) turns into (D — 4). 

Finally, terms involving get a factor J2 n =i Z in ^nO = -/>A (4) owing to JZZS) and (STEty . 
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Inserting the Lorentz decompositions of the tensor integrals into ()6.10|) . we find the 
following reduction equations for the tensor coefficients upon comparing coefficients of 
covariants, 
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Since we have distinguished the index k in the derivation of ()6.12|) . the resulting tensor 
coefficients P... are not symmetric under the exchange of k with one of the indices i r , 
r = 1, . . . , P. In order to distinguish them from the symmetric tensor coefficients E___, we 
marked them with a bar. Symmetric tensor coefficients can be easily obtained by adding 
all P results with k exchanged with one of the i r and dividing the sum by P, e.g. , 
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In ()6.15|) . E 00il and E 0il0 are determined from ()6.13|) . while E il00 is determined from 
(EH- 

For P > 4 extra terms of order (D — 4)P o... have to be added to the equations ()6.12|) 
and ()6.13|) . For P = 4 the last-but-one contribution in ()6.1()|) is of 0(D — 4), but the last 
term yields a finite contribution for D — > 4, because the coefficient Poooooo is UV divergent. 
We calculate this contribution upon inserting £ ,a/3Ml ' M4 |div = {ggg}" 13 ^ 1 "^Pooooooldiv into 
(|6.4|) and using (D — 4)£ , ooooo from (|A.5|) . After symmetrizing in the Lorentz indices, we 
get 
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This contribution to the coefficients -Ejoooo can be included by replacing — D 0000 (0) in 
(16.12}) by — XjQ Doooo(O) + ^ • The cases P > 4 can be treated analogously, but usually 
do not appear in renormalizable quantum field theories. 

After the symmetrization, we thus find for the tensor coefficients up to rank 5: 
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For the 4-point tensor coefficients that result from omitting iVo in the 5-point integrals, 
we have introduced the auxiliary quantities 
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and similar quantities resulting from these relations with index pairs "00" added to the 
A..(0) functions on both sides. 
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7 Reduction of 6-point integrals 



Following the guideline of the reduction of the scalar 6-point integral to six scalar 
5-point integrals [123, the 6-point tensor integrals of rank P can be reduced to six 5-point 
tensor integrals of rank P as described in Ref. [ |^] . This method, which was used in the 
calculation of one- loop corrections to e + e~ — > 4/ [EE], is more explicitly worked out in 
App.EJ 

In the following we describe a method that reduces 6-point tensor integrals of rank P 
to 5-point tensor integrals of rank (P — 1). The scalar 6-point integral should be treated 
following Refs. [I22ES] as explicitly described in App. [0J The tensor reduction can be 
derived by considering the determinant 
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(7.1) 



The r.h.s. is obtained by adding the (k + l)th row to the first row and using ([5.4)1 . 

In four dimensions, this determinant vanishes, as can be seen from the first form in 
(JUIJ), because q is linearly dependent on the four (non-exceptional) momenta pi, i = 
1,...,5, % 7^ k. We again do not use this fact, but translate the integral over T into 
a form that has a factor of 0(D — 4) rendering the whole contribution zero for finite 
integrals. Inserting (|7.1|) into the integrand of the tensor integral F fll '- fJ-p results in 
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We expand the determinant along the (k + l)th row and use the fact that the four- 
dimensional metric tensor can be written as 
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for two arbitrary sets of linear independent momenta Pi,P2>P3>P4 an d fci, /c2, ^3, &4- This 
yields 

'r = -Xj$Fw..*rQ£- g(A *). (7 . 4) 



Inserting the r.h.s. of (|7.1|) into the integrand of the tensor integral i^-Mp results in 



2pipi 



A' 
Pi 



Pjfc-1 



Pk+1 



P5 



2pfe-ipi 
/1 

2pfc+iPi 
2p 5 Pi 



2pip 5 

2pfc-lP5 

/ 5 

2pfc+lP5 

2p 5 P5 



(7.5) 



Expanding the determinant along the first row and the first column according to the 
analogue of (|2.22|) . yields 



J* 7 



(fcm)(0n)^m 

n,m=l 



i-.-mp 



(0) 



(7.6) 



where ^(Lvc M is S iven in (EH- 
From (|7.4|) and ()7.6|) we obtain 



5 

A (H(o«)^ m 

n,m=l 



l— MP 1 



71 



^•^(0) +xf ) F^-^(^-^ (4) ^ 



(7.7) 

The last term in ()7.7|) only contributes in dimensional regularization if F^^-vp is singular. 
For UV singularities this is the case if P > 7, which is usually not needed in renormalizable 
theories. As explained in Section IB~Hl IR (soft and collinear) singularities of F afJ-1 '" fJ-p only 
appear in contributions that are proportional to a momentum pf. These contributions 
vanish exactly in (|7.7jl . Therefore, the terms involving (g — g^) in (|7.7J1 can be omitted 
for P < 7. 8 For P > 7 the inverse determinant that is implicitly contained in gu) can 
always be cancelled. 9 



8 This is again in agreement with the observation [ 29 that in the absence of UV divergences reduction 
formulas valid in 4 dimensions remain valid in D dimensions up to terms of 0{D — 4), independent of 
possible IR singularities. 



9 According to 



X 



(5) 

fed 



YTn=i Zfcnfn- F° r each of these terms, Z^gu) can be expressed via 



(5), 



(|7.3() by a determinant without denominator. 
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Introducing the matrix 



M, 



(k) 



( 2piPi 

Ipk-iPi 
fi 

^Pk+iPi 

\ 2P5P1 



%PlP5 \ 
2Pk-lP5 

h 

^Pk+Wb 

2P5P5 J 



(7- 



(17. 7|) can be written as 



^£tl...£tp 



E E (M7k) 



71=1 m =! 



fm 



E 



n) 



l—HPI 



(o) +F a ^(gZ-g(^ y ), (7.9) 



which expresses the 6-point tensor integral of rank P in terms of six 5-point tensor integrals 
of rank (P — 1). The inverse of Mn* is given by 



K) 1 



-X 



(5) 



(5) 



(fcj)(0i)/^ v «) , 



i,j,fc = 1,...,JV. 



(7-10) 



In the form (|7.9|) our result can easily be extended to the reduction of iV-point functions 
with iV > 6 by simply forming a matrix similar to Mha by selecting five momenta for 
the columns and four momenta for the rows out of the (N — 1) available momenta of the 
iV-point function. 

Equation (|7.7|) can also be used to derive an alternative reduction of tensor 6-point 
integrals. Multiplying it with Xko, summing over k = 1,...,N, and using f)2.30|) . for 



X 



(5) 
(km)(0n) 



yields 



det(X (5) )F m = E X nh ) m\ E '" 1 "' IXp { n ) -E^-^iO) 



n,m=l 



+ det(X®)F a ^(gZ-g { 4 



(7.H) 



Here, as in (j6.10|) . all inverse Gram determinants have been absorbed in the four-dimensi- 
onal metric tensor, which appears only in the difference (g — gu))- The result (|7.11j) is 
equivalent to Eq. (64) of Ref. [H2J. 

Finally, we insert the decompositions of tensor 6-point integrals into Lorentz covari- 
ants in order to derive explicit reduction formulas for the tensor coefficients. Since we 
consider only tensors up to rank 3, we can omit the terms involving (g — gu\). The tensor 
decompositions explicitly read 

^ = E p£*ix> F " u = E + 9" U F 00 , 

U=l 11,22=1 



f» vp = E ptiPl^Fm^ + EWV 

*l,«2,i3 = l H=l 



(7.12) 
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In four dimensions, some covariants in these decompositions are redundant in the sense 
that they can be expressed by the others. For instance, in the decomposition of F M one of 
the five covariants p^F^ is redundant, because one of the momenta p^ can be expressed 
by the other four linearly independent vectors. Similarly, all covariants involving metric 
tensors are redundant. However, by keeping these coefficients we can avoid the appearance 
of explicit inverse Gram determinants in the reduction formulas. 

Inserting the Lorentz decompositions of the tensor integrals in the reduction formu- 
las given above, we can read off the reduction formulas for the tensor coefficients upon 
comparing coefficients of covariants on both sides. Generically we find 



p.. . 
- 1 ji\...ip 



5 

n=l 



E (il)u-(ip)n( n ) S hn ■ ■ -bipn ~ E h ... ip (0) 



P < 7, 



(7.13) 



with 

Con — Cfcn = 0, Cj n = (jtfy^j , j, 

for the reduction given in (|7.9jl and with 



n 



. , 5, j ^ k 



c 0n = 0, c jn = xgV det(X( 5 )) = (x< 5 >) l n , j, 



n 



1,...,5 



(7.14) 



(7.15) 



for the reduction given in ([7.110 . In the numerical reduction we can select the equation 
that is numerically most stable. For example, in ([7.14[l we can choose k such that the 

~ (5) 

modulus of Xj^Q = — det M( fc ) is maximal. 

Since we have distinguished one momentum in the derivation of (17. 9[) the resulting 
tensor coefficients F are not symmetric under the exchange of j with one of the indices 
i r . This can be easily cured as in the case of 5-point functions [see ([6.14(1 and ([6.15(1 ] by 
adding all P results with j exchanged with one of the i r , and dividing the sum by P. 

Thus, we find from ([7.13(1 for the tensor coefficients up to rank 3 



U«2 



OOii 



<k in [E (n)-E (0) 

n=l 

o, 

1 5 



% x = 1, . . . ,5, 



n=l 
5 



- J2\ c hn E {i2)n (n)5 i2n - E i2 (0) 



J2 dm E 00 (n) - E 00 (0) 



n=l 
5 



n / y ) c im E(i 2 ) n (i 3 ) n (n)5i 2n 5i 3n -£/j 2 j 3 (0) 

6 n=l L 

ix,i2,k = 1, ... ,5. 



+ (h «-> H) + (h «-> H t 



(7.16) 



+ (*i <-> ii) k ii,i2 = 1, • • • ,5, (7.17) 



(7.18) 



For the 5-point tensor coefficients that result from omitting Nq in the 6-point integrals, 
we have again used the auxiliary quantities 



^ 1 (0)=£; il _ 1 (0), ii = 2,...,5, 
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£i(0) 



-$> n (0)-£ (0), 



(7.19) 



n=2 



#U-l,i 2 -l(0), *1j *2 = 2, . . . ,5, 



5 



£ln(0) 



-]T£ mi (0)-£ n (0), i 1 = l,. 



,5, 



(7.20) 



n=2 




^ii-l.ia-Ms-lCO)) ^1,^2,^3 — 2,. 



,5, 



5 



-E , liii 2 (0) 



- ^Knj 2 (0) - -Eni 2 (0), h,i 2 



1,...,5. 



(7.21) 



n=2 



8 Summary 

Methods for a systematic evaluation of one-loop tensor integrals have been described 
for graphs with up to six external legs. The results are presented in a form that can be 
directly translated into a computer code; only the scalar 3- and 4-point integrals have to 
be taken from elsewhere. 

While UV divergences are treated in dimensional regularization, possible IR (soft or 
collinear) divergences can be regularized either dimensionally or with small mass param- 
eters; the described results are valid in either IR regularization scheme. Moreover, the 
results hold if internal masses are complex parameters, which naturally appear for unsta- 
ble internal particles. The generalization of the proposed methods to functions with more 
than six external lines is straightforward. 

Particular attention is paid to the issue of numerical stability. For 1- and 2-point inte- 
grals of arbitrary tensor rank, general numerically stable results are presented. For 3- and 
4-point tensor integrals, serious numerical instabilities are known to arise in the frequently 
used Passarino-Veltman reduction if Gram determinants built of external momenta be- 
come small. For these cases we have developed dedicated reduction techniques. One of 
the techniques replaces the standard scalar integral by a specific tensor coefficient that 
can be safely evaluated numerically and reduces the remaining tensor coefficients as well 
as the standard scalar integral to the new set of basis integrals. In this scheme no danger- 
ous inverse Gram determinants occur, but inverse modified Cayley determinants instead. 
In a second class of techniques we keep the basis set of standard scalar integrals and 
iteratively deduce the tensor coefficients up to terms that are systematically suppressed 
by small Gram determinants or by other kinematical determinants in specific kinematical 
configurations. The convergence of the iteration can be systematically improved upon 
including higher tensor ranks. For 5- and 6-point tensor integrals, we describe reductions 
to 5- and 4-point integrals, respectively, that do not involve inverse Gram determinants ei- 
ther. Compared to some other existing methods, the described methods are distinguished 
by the fact that the reduction from 6-(5-) to 5-(4-)point integrals decreases the tensor 
rank at the same time. 

We finally emphasize that the presented methods have already been successfully ap- 
plied in the calculation of a complete one-loop correction to a 2 — > 4 scattering reaction, 
viz. the electroweak corrections to the charged-current processes e + e~ — > Af. The de- 
scribed methods, thus, have proven their reliability in practice and will certainly be used 
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in future loop calculations for interesting many-particle production processes at the LHC 
and ILC. 
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Appendix 

A UV-divergent parts of tensor integrals 



In the reduction formulas given above, products of (D — 4) with tensor integrals 
appear. These give rise to finite terms originating from UV singularities in the loop 
integrals. As mentioned above, no IR-singular integrals multiplied with (D — 4) appear in 
the reduction formulas. The UV-singular parts of the loop integrals can be derived easily 
from the Feynman-parameter representation or by using (|5.1Uj) for these parts only. In 
the following, we list results for (D — 4) times one-loop integrals omitting terms of order 
0(D — 4). For the 1-point functions A ..(mo) we get 

m 2n+2 

(D - A) A = -2ml ( g - 4 )^ = - 2 n-i(°n + l)! ' " = 1,2 ' ' ' ' (A ' 1} 

2n 

For the IR-finite 2-point functions B.„(pi, mo, m i), i-e. excluding the case p\ = m\ = 
ml = 0, we obtain 



(D - 4) B 
(D - 4) B 1 
(D - 4) B 00 
(D - 4) £?ooi 
(D - 4) Eoooo 
(D - 4) Soon 
(D - 4) Sooooi 
(D - 4) Boom 



-2, 



|(p 2 -3m2-3m 2 ), (D — 4) B n = - 
i_(p 2 -2m 2 -4m 2 ), (D - A) B ni 



2 

3' 



12 

1 

120 



1 

2' 



60 

1 
240 



Pi ~ ^Pi\ m o + m i) + 10(m + m Q m 1 + m-J 
15m 2 ), (D-A)B nn = 



5ml 



2 

5' 



p\ — Ap\m\ — 6pfmf + 5mQ + lOm^mf + 15m 



,2_2 



2_2 



1 
3" 



(A.2) 



-±{2p\ - 3m 2 - 12m 2 ), (D - 4) B nm 
For the 3-point functions C...(pi, p 2 , m , m l5 m 2 ) we obtain, denoting (p x — p 2 ) 2 = s 12 , 
(D - 4) Coo = 



l 

' 2 : 
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(D - 4) C 00l 

(L>-4)C 000 o 
(D - 4) C 00ii 



1 

6' 



^[Sl2 + p\ + Pi] ~ J2( m l + m l + m 2) ? 

(£> — 4) Cooij — — ^, 



j_ 

12' 



(£>-4)C, 



OOOOi 



1 

'240 



2s 12 - 5m^ + _ 5m 2)(! + ( j. n ) 

n=l 



(£> — 4) C om — ^j, (D — 4) Cooiij — ^, 



(D-4)C ( 



oooooo — 



2880 



2s\ 2 - Qs 12 m 2 + 30m^ + 2s 12 ^(Pn _ 6m n) 

n=l 

6r^E(2^-5m^) 

n=l 

2 

Z (PmPn - 6p^m^ + 15m 2 m m 2 n )(l + 5 mn )] , 



m,n=l 



1 

720 



(-D — 4) Cooooii 

(D — 4) Cooooij = 7I0 
(-D — 4) Coouu 



3s 12 - 6m^ + ^(p^ _ 6m 2 )(1 + 25m) 

n=l 
2 

2si 2 - 3mg + (Pi - 6m l 



n=l 



J_ 

30' 



(D - 4) Coouij — — lio' ~ 4) Cooiiji — _ t|o' 

2 



(D - 4) Cooooooi = YoOso 3s i2 - 7si 2 mQ + 2\m\ + s i2 Z _ 7m «)( 2 + <*»n) 



n=l 



(D — 4) CoooOiii — ~TfiS0 



(D — 4) CoooOiij 



-7mgXK-37r£)(l + «f in ) 

n=l 

2 

+ Z (PmPl ~ 7 Pm m l + 21m 2 m m 2 n ){l + 25 im 5 in ) 
=i 

4si2 _7 m 2 + ^ (p 2_ 7m 2 )(l + 35 . n) 



m,n=l 
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n=l 
2 



6si2 - 7m 2 + J2(pI- 7m l)( 2 + <Ji„) 



n=l 



(D — 4) Cooniii — (D — 4) Coouuj — 2T0' (-^ ~~ ^) Coomjj — 4^5 > (A-3) 

where i,j = 1,2 but i 7^ j. All other 3-point tensor coefficients up to rank 7 are UV 
finite, so that for them (D — 4)C... = if they are IR finite. For the 4-point functions 
D...(p 1 ,p 2 ,p 3 ,m ,m 1 ,m 2 ,m 3 ) we find, denoting (p l - p 2 ) 2 = s 12 , (pi - p 3 ) 2 = s 13 , and 
(P2 - P3) 2 = s 23 : 



(D — 4)D 



0000 
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- 4) -Doooooo = 


(D- 


— 4) -DoOOOii 


(D- 


' -DooOOOOi = 


(D- 


" 4) -Dooooin — 
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i|o[ s 12 + Sl3 + S 23 + p\ + p\ + Pi] ~ M m O + m l + m 2 + m s)> 

_ i2o> (- ~~ 4) Dooooij = -2M), 



2880 



Pn{ 1 + din) + S mn (l + S in + 4 



.n=l m,n=l 
m>n 



^ 480 

n=0 



2io> (D — 4) -Doooomj — 7^' (-^ — 4) D'ooooijfc — j^q, (A-4) 

where z, j, fc = 1, 2, 3 but are pairwise different. All other 4-point tensor coefficients up to 
rank 7 are UV finite. 

For the 5-point functions E.„(pi,P2,P3,P4, m , mi, m 2 , m 3 , m^), there is only one UV- 
singular tensor coefficient up to rank 6, 

(D- 4) Doooooo = (A.5) 

B Tensor coefRcients of singular 3-point functions 

The vanishing of the modified Cayley determinant det(X^), as defined via ()2.25|) . 
is a necessary condition for the existence of a leading Landau singularity in a one-loop 
TV-point integral. For 3-point integrals this means that det = for IR-singular (either 
soft or collinear) integrals, so that the reduction methods of Sections 15.21 and 15.31 are not 
applicable in this case. If in addition the Gram determinant is small, for IR-singular 
3-point integrals also the X j are small, and the reduction method of Section 15.41 cannot 
be used either. One could still, however, use the method of Section 15.51 In the following 
we describe a way of evaluating these specific 3-point functions that does not make use 
of an iteration technique, but is based on analytical simplifications that are admitted by 
the simple structure of the special cases. 

The simplifications are achieved by directly using the analytical results for the standard 
scalar integrals and for the tensor coefficients, as obtained with the Passarino-Veltman 
reduction, and by rewriting them in such a way that the limit of vanishing Gram de- 
terminant does not involve numerical cancellations. To this end, the scalar integrals are 
split into two parts: one contains the asymptotic behaviour of the integral in the limit of 
vanishing Gram determinant A® up to a specific order n and a corresponding remainder 
which is of O (JA( 3 )]™ +1 ^ . We symbolize this splitting by introducing the asymptotic op- 
erators T^} and TVj£l X0 , which define the asymptotic behaviour of a function f(x) for 
x -> x by 

/(*) = t£Lo [/(*)] + [/(*)] > [/(*)] = ° - *°) n+1 ) > n = °> • • • ■ 

(B.l) 

If the function f(x) is analytical at x = xq, X^} Xq is the usual operator for a Taylor 
expansion up to order n. 

Making use of these definitions, we now describe the treatment of the IR-singular 3- 
point tensor integrals that were needed in the calculation of the one-loop corrections to 
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e + e — > 4/ It is convenient to switch from the original definition (|2.1|) of arguments 
on tensor coefficients to the new notation 

B...(pl,m ,mi) = B„.(p 1 ,m ,mi) ) 
C...(pl, (pa - Pi) 2 ,pl,m ,m 1 ,m 2 ) = C..,(p 1 ,p 2 , m , mi, m 2 ). (B.2) 

(%) Collinear- singular case with two off-shell legs: C,..(m 2 , s, s', 0, m, M) 

Here m denotes a small real mass, which will be neglected whenever possible. In this 
limit the relevant scalar integrals read 



B (0) = B (s, m, M) = A + In(^) + 2 + {— - l) ln( 

+ 2 + ( — - l\ In 
B (2) = 5 (m 2 , 0, m) = A + ln(-Q + 2, 
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In 



2 

Vm 2 , 
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M 2 - s' 
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M 2 — 

M 2 



\M 2 - s' 



7 +Li 2 



^-LiJ^ 



M 2 



VM 2 



(B-3) 



where M 2 is complex with a finite or infinitesimal negative imaginary part, which is also 
present for vanishing M 2 . The Gram determinant is given by 



A (2) 



J\2 



(B.4) 



so that the delicate limit is 5s = s' — s — > 0. The asymptotic expansions of the scalar 
integrals in ()B.3j) for this limit can be worked out easily; the first few terms read 
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5s 
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s-M 2 



M 2 , 
1 + ln| 
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5s 



(^ 1 )] +K ^»i B »< 1 ») 



+ 2- 
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2s(M 2 - s) 
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+ Kfo-to [Co 



/M 2 -s\ M 2 
V m 2 J s 



5s(M 2 - 2s) 



2s(M 2 



In 



M 2 



M 2 



(B.5) 



where we have kept s fixed. Inserting these or forms with more explicit terms of the 
asymptotic expansion for the scalar integrals into the explicit formulas for the tensor 
coefficients, one obtains expressions like 

_ M 2 [(M 2 - s)s - 5s(AM 2 - 5s)} M 2 - s + 5s ( M 2 — s 



2s 2 (M 2 - s) 2 



2(M 2 - s) 2 



M 2 [M 2 s{M 2 -s)- 5s(4M 4 - 7M 2 s + 2s 2 )} M 2 - s 



2(M 2 - s) 2 s 3 



M 2 



2{s + 5s) 
{5s) 2 



n 



(2) 

Ss^O 



[Bo(l)} 



M 2 



s — 5s 



5s 



-n 



(i) 



[Co], 
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M 2 J s 5s' 



1. 1, / /i 2 \ M 2 (M 2 — 5s) (M 2 -s\ M 2 + 2s-5s 

_M ! ^ £ Bo(i)] 
4os 

Here the orders n in the 7?A n ) operators are chosen in such a way that all terms involv- 
ing IZ^ contribute only in 0(8s) in spite of the enhancement factors 1/Ss m . Note that 
no delicate cancellations for 5s — > appear in the other terms, although the original 
Passarino-Veltman results contain plenty of terms involving l/5s m in front of linear com- 
binations of scalar integrals. Thus, the above forms are numerically stable as long as 
the remainder terms 7lS n ) can be evaluated in a stable way. This task is, however, eas- 
ily achieved upon expanding the scalar integrals as in (|B.5J) to a high order, e.g., with 
computer-algebraic methods, and dropping the first n orders. The resulting series are easy 
to evaluate, and an arbitrarily high precision can be achieved by including sufficiently high 
orders in the expansions. On the other hand, if 5s is not small, the TZ^ terms can safely 
be evaluated upon numerically subtracting the terms from the scalar integrals. In 
this way an arbitrarily high precision can be achieved as long as s, s' ^ and s ^ M 2 . 
The case s = M 2 does not occur in our application, the cases s = and s' = are treated 
below. 

(ii) Collinear- singular case with one off-shell leg: C...(m 2 , 0, s', 0, m, M) 
Specializing the previous case to s = 0, the scalar integrals read 

50(0) = B (0, m, M) = A + ^(jp) + ^ 

«- L M£H^) (B - 7) 

with B (l) and B (2) still as given in (|B.3|) . The limit of vanishing Gram determinant is 
now reached for s' — > 0, where the scalar integrals can be expanded according to 



s> 



5 (l) = 5o(0) + ^ + 7^lo [B (l)], 



c = 1 



M 2 



s' \ , fM 2 \ , s' 1 <d 
l + ^-=7 ln(— =- +1 + ^ () 



2M 2 J \m 2 J AM 2 



+ KLo [Co] , (B.8) 



or to higher orders if needed. Making use of these expansions, the first few tensor coeffi- 
cients can be written as 



X . 1 , / u 2 \ 3M 2 + s' M 2 — s' (1) .„ . . . 
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The IZ^ terms, which are suppressed by a factor (s') n+1 , can be evaluated to arbitrary 
precision for all values of s' as described above. 

(Hi) Collinear- singular case with one off-shell leg: C,,,{m 2 , s, 0, 0, m, M) 
Specializing case [(I)] to s' = 0, the scalar integrals read 

t 
M 2 



So(1) = S (0,0,M) = A + 1b(^ ? )+1, 



with -Bo(O) and -Bo (2) still as given in ()B.3j) . The limit of vanishing Gram determinant is 
reached for s — > 0, where the scalar integrals can be expanded according to 



1 r/l+ 2^)K^) +1 -i| J ] + ^^- 

or to higher orders if needed. Making use of these expansions, the first few tensor coeffi- 
cients can be written as 

° 2 = ^p + Po(0)] , 



M 2 



C »=4 A+ 4n^J + ^ 5 -'C[*(0)]. (B.12) 

The TZ^ terms, which are suppressed by a factor s n+1 , can be evaluated to arbitrary 
precision for all values of s as described above. 

(iv) Soft-singular case: C...(mf , s, m 2 ,, A, mi, m 2 ) 

For processes with external fermions in the massless limit (rrii — > 0), the Passarino- 
Veltman reduction of this case turns out to be less delicate than the previous ones. In fact, 
no special treatment was necessary for e + e~ — > 4/ [ESI, although one could also improve 
the stability as described in the previous sections. We attribute the robustness of this case 
to the following reasons. Firstly, because A is an infinitesimal photon mass fx = fi = 
and all 3-point tensor coefficients are directly obtained from 2-point coefficients without 
further recursions. Thus, instabilities do not accumulate. Secondly, for massless fermions 
the Gram determinant A^ 2 ^ = — s 2 vanishes only for s — > 0, and this case appears for 
e + e~ — > 4/ only in regions of phase space that are suppressed by r w /M w . 

C Alternative reduction of 5-point integrals 

In Ref. [ 1213 we have worked out a reduction of 5-point tensor integrals that follows the 
strategy proposed in Ref. [122] for scalar integrals in four space-time dimensions. Here 
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we briefly describe the derivation of this method in D dimensions, to make closer contact 
to the methods used in this paper. 

The reduction is based on different ways of evaluating the determinant 



2q 2 2qp 1 
2p x q 2pipi 

2p A q 2p A p x 



2qp A 

2pm 

2p A p A 



(C.l) 



which vanishes in four dimensions owing to the linear dependence of any five momenta. 
In D dimensions the integral over £' can be easily evaluated to 
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A D Q 
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2g a f3 2p ljQ 

2p x ,p 2p x p x 



2p A)l3 2p A p 1 ... 



2p4,a 

2p A p A 
2p A p A 



(C.2) 



where we have identified the form ()6.3|) of the metric tensor gu), a p m f° ur dimensions. 

On the other hand, the integral over £' can be evaluated in terms of 4-point functions 
as described in Section 2 of Ref. [|2S| with the only difference that no additional UV 
regularization is needed, because we now keep the dimension D general. In detail, this 
means that the factor —A 2 /(q 2 — A 2 ) introduced in (2.5) of Ref. is absent, and the 
result analogous to (2.19) of Ref. [ 26J becomes 



det(Y)£ w " 



HP 



J2det(Y n )D^(n)+ £ Z^2p m , a V a ^^ 



n 



n=0 



n,m=l 



(C.3) 



where Y = (Yij), i,j = 0, . . .4, was defined in (|2.26jl . and Y n is obtained from the 5- 
dimensional modified Cayley matrix Y by replacing all entries in the nth column by 1. 
The last term of (j(J.3|) . which results from (jC2|) . contributes only if ^pm-mf involves 
a divergent coefficient -Eoo... corresponding to a covariant containing a metric tensor. As 
explained in Section 15.81 such coefficients are free of IR divergences, and power counting 
shows that UV divergences only occur for P > 4. Therefore, the last term in (|C3|) is of 
0(D — 4), and thus irrelevant, for P < 3. For P = 4, this term can be explicitly evaluated 
using (jA.5|) yielding 



det (Y)E' 



Ml- --£14 



J2det(Y n )D^(n)+ £ Z^p^V^" 

n=0 n,m=l 

- — A (4) W" 1 -" 4 —J 

48 iyy; 24(£> - 4) 



■fa i 



n 



A^{(g-g {A) )gr^\ 



(C.4) 
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where {(g — 5 , (4))g} Atl is a symmetric tensor of rank 4 constructed according to the 
rules explained in Sectional 



{(g-9(4))g} 



" , " IH = (5 - g^Y^g^ 4 + (g- g^Y^g^ + (g- g^Y^g^ 3 
+ {g- g(i)Y 2tls g^ 4 + {g- g^)Y 2tH g^ 3 + (g- g^Y^g^ 2 - 



(C.5) 



The first term in the last line of (jC4|) is just the finite contribution defined in 

(2.15) of Ref. an d the UV-divergent terms of the second term in the last line exactly 
cancel the UV divergences of the 4-point integrals in the first line. Thus, the result for 
£W-M4 exactly receives the form of (2.19) of Ref. [EE!, 



det(F)^ 1 -^ = -Y, det{Y n ) D^^-^in) + £ 2p m ^V^ a ^-^{n) 

n,m=l 



n=0 



(C.6) 



where the superscript "(fin)" indicates that the UV parts have to be consistently omitted, 
as e.g. following the MS prescription. 



D Alternative reduction of 6-point integrals 

Here we describe the reduction of 6-point tensor integrals of rank P (including the 
scalar case P = 0) to six 5-point tensor integrals of equal rank that is based on the strategy 
of Ref. [123- This reduction is related to the reduction of 5-point functions as given in 
Ref. [EHj and App. [njand has been used in the calculation of the electroweak corrections 
to e + e + — > 4f [EH)- Moreover, it is needed to reduce the scalar 6-point function to 5-point 
functions [E2j- 

It starts from the observation that 



4-D 



17T 



J iVoiVi • • • 



7 mp 



N + Y 00 2q Pl 
Yio - loo 2pipi 



2pip 5 



Ygo - Y 00 2p 5 pt . . . 2p 5 p 5 



(D.l) 



which is correct in any space-time dimension D as long as the five four-momenta pi 
(i = 1, . . . , 5) are linearly dependent, and thus for four-dimensional pi, because then the 
five last columns of the determinant are linearly dependent for an arbitrary D-dimensional 
momentum q. The l.h.s. of this relation is practically the same as in Eq. (2.10) of Ref. [I2"S]. 
where the reduction of 5-point integrals is described. The same manipulations as described 
there lead to the result 



pui—np 
1 
1 



-£w-w(0) 

loo 
Y w 



>oi 



Yt 
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Y 
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Y05 
Y 15 



Y f 
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Equation (jD.2|) expresses _F M1, " Mp in terms of six 5-point integrals, 

F h..*p = Vn E^(n) with ^ = -±21, (D.3) 

n=o aet{Y ) 

where Y = (ly), i, j = 0, . . . 5, and Y n is obtained from the 6-dimensional modified Cayley 
matrix Y by replacing all entries in the nth column by 1. For the scalar integral Fq, this 
result is identical with the one of Ref. [ 22J. 

By inserting the Lorentz decompositions as given in ([7.12)1 . we can derive explicit 
formulas for the scalar 6-point function and the coefficients of tensor 6-point integrals 
from flUU): 

5 

F = -$> n £ (n), (D.4) 

n=0 
5 

Fh = - VnE {il)n (n)5 nn - 770^(0), *i = 1, - . . , 5, (D.5) 

n=l 
5 

Foo = VnE 00 (n), 

n=0 
5 

F hi 2 = -Y,^E{i 1 ) n ^ 2 ) n (n)6 iin 6 i2n -r]oE ili2 (0), i u i 2 = 1, . . . , 5, (D.6) 

n=l 
5 

Fooh = -J2VnEo (i 1 ) n {n)5 hn -i] E oi 1 {0), i x = 1, . . . , 5, 

n=l 
5 

-^1*2*3 = _ 7 7r l -E , (i 1 )„(i 2 )„(i3)„( n )^in^ 2 n5i 3 n ~ T)oEi li2 i 3 (0) , «1, «2, «3 = 1, • - • , 5. (D.7) 
n=l 

The 5-point tensor coefficients that result from omitting iVo in the 6-point integrals have 
been given in (j7.19J) . 
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